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Abstract 

We construct a perturbative S-matrix for interacting massive scalar fields in global 
de Sitter space. Our S-matrix is formulated in terms of asymptotic particle states in 
the far past and future, taking appropriate care for light fields whose wavefunctions 
decay only very slowly near the de Sitter conformal boundaries. An alternative for- 
mulation expresses this S-matrix in terms of residues of poles in analytically-continued 
Euclidean correlators (computed in perturbation theory), making it clear that the 
standard Minkowski-space result is obtained in the flat-space limit. Our S-matrix 
transforms properly under CPT, is invariant under the de Sitter isometries and pertur- 
bative field redefinitions, and is unitary. This unitarity implies a de Sitter version of 
the optical theorem. We explicitly verify these properties to second order in the cou- 
pling for a general cubic interaction, including both tree- and loop-level contributions. 
Contrary to other statements in the literature, we find that a particle of any positive 
mass may decay at tree level to any number of particles, each of arbitrary positive 
masses. In particular, even very light fields (in the complementary series of de Sitter 
representations) are not protected from tree-level decays. 
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1 Introduction 



The S-matrix is an invaluable tool for studying quantum field theory (QFT) on Minkowski 
space. Even leaving aside the all-important connection to experiments, its gauge invariance 
and invariance under field redefinitions make the S-matrix a powerful way to organize our 
understanding of quantum fields. The Coleman-Mandula [1], Haag-Lopuszanski-Sohnius 
[2], and Weinberg- Witten [3] theorems are prime examples of the utility of an S-matrix 
approach to field theory, and the above properties allow the S-matrix to be well-defined even 
in perturbative (and perhaps also nonperturbative) string theory. 

At a more mundane level, the fact that all correct calculations of an S-matrix element 
must agree allows a clean comparison of different approaches, choices of gauge, etc., that 
helps to resolve potential controversies and hastens the advance of knowledge. In contrast, 
the lack of an S-matrix-like object has been sorely felt in years of controversy regarding 
quantum fields in de Sitter space. Arguments continue over the interpretation of calculations 
in both simple self-interacting scalar theories on a fixed de Sitter background (compare e.g. 
[4, 5, 6, 7, 8, 9, 10]) and also in the more complicated case of gravitational theories (compare 
e.g. [11, 12, 13, 14, 15, 16, 17]). 

In this work, we introduce an S-matrix for weakly-coupled quantum field theories in de 
Sitter space that can be computed order-by-order in perturbation theory. The associated 
spaces of asymptotic states are defined using the interacting Hartle-Hawking state as the 
vacuum. For theories of massive (M 2 > 0) scalar fields, we show our S-matrix to be unitary, 
de Sitter invariant, and invariant under perturbative field redefinitions. It also transforms 
properly under CPT and reduces to the usual S-matrix in the flat-space limit. Our analysis 
is strictly perturbative and we consider in detail only theories of interacting scalars. Our 
final discussion (section 8) will comment briefly on extensions to gauge fields and why gauge- 
invariance is to be expected. Perhaps our construction will also be of use in perturbative 
string theory on (likely meta-stable) de Sitter backgrounds. 

To be specific, in this work we use the term a de Sitter (dS)" to refer to global de Sitter 
space, including both the contracting and expanding cosmological regions; see figure 1. Thus 
we formulate a global de Sitter S-matrix below, though we believe that an analogous S-matrix 
can be defined for both the Poincare patch (also known as the k = cosmological patch), 
and the hyperbolic (A; = — 1) cosmological patch of dS. Furthermore, these S-matrices should 
all be closely related through appropriate analytic continuations 1 . 

One sometimes hears the claim that there can be no S-matrix on de Sitter space. We are 
aware of the following concerns regarding potential S-matrices: 

i) The Minkowski S-matrix is defined using in/out perturbation theory, but it is well- 
known that in/out perturbation theory in dS suffers from infrared (IR) divergences. So 
this definition does not work in de Sitter space. 

ii) There is no positive-definite energy-like conserved quantity in de Sitter space. As a 

1 While there should also be an S-matrix for the static patch, it would naturally be denned using asymptotic 
particle states built on the static vacuum in contrast to the Hartle-Hawking state used in our construction. 
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Figure 1: Left: A finite piece of global D-dimensional de Sitter space represented as a 
timelike hyperboloid in D + 1 Minkowski space. The diagonal line is a cosmological horizon. 
Center: A conformal (Carter-Penrose) diagram for global de Sitter marked with horizontal 
cross-sections, each representing an The left and right edges are the poles. Right: 

The dS conformal diagram showing a cosmological horizon %. Each point on the diagram 
represents an S D ~ 2 which contracts to zero size at the left and right edges. 



result, one-particle states can decay and all particles are unstable; thus there are no 
viable asymptotic states. This is directly related to the concerns of e.g. [18, 19, 20, 21]. 

iii) The contracting phase of global de Sitter space tends to blueshift particles to high 
energies. In a theory with dynamical gravity, many states which are weakly-coupled near 
past infinity induce large gravitational back-reaction near the minimal-radius sphere (the 
de Sitter "neck," rj = in the coordinates of (2.2) below). Semi-classically, this should 
result in gravitational collapse to a cosmological singularity. There is then no reason 
to expect late-time behavior described by weakly-coupled asymptotic states near the 
future de Sitter boundary. 

iv) At least in string theory, all known de Sitter vacua are at best met a-st able. So one ex- 
pects that mere particle excitations of a de Sitter background cannot provide a complete 
set of outgoing states. 

v) The causal structure of global de Sitter space, and in particular the fact that its 
past/future boundary is spacelike, prevents any one observer from interacting with a 
complete set of ingoing/outgoing states. This means that the S-matrix is not experi- 
mentally accessible to a single observer and hence need not necessarily be a well-defined 
object in a fundamental theory. 

The reader will note that comments (ii)-(v) are not directly relevant at the level we wish 
to work. This is particularly manifest for issues (iii), (iv), and (v), which concern dynamical 
gravity, string theory, or other supposed fundamental theories. Issues (iii) and (iv) (and 
arguably also (v)) are also intrinsically non-perturbative, and so do not obstruct the more 
modest goal of formulating a perturbative S-matrix even for gravity or string theory. While 
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issue (v) implies that a direct connection to experiment is unlikely, a de Sitter S-matrix 
should nevertheless provide a useful theoretical tool along the lines noted above. 

Let us briefly elaborate on the irrelevance of issue (ii) , the expected decay of all particles 
in de Sitter space. First recall the usual situation in flat space: at the perturbative level, 
each free field is associated with an appropriate set of asymptotic particle states, whether or 
not these particles turn out to be stable. In the unstable case there is a nonzero l-to-2 (or 
more) amplitude that can be computed order by order, and this amplitude is related by the 
optical theorem to the self-energy correction for that particle. While this implies certain IR 
divergences in the computation of the strict order-by-order S-matrix, 2 the physics of these 
divergences is well-understood. We will show that all of this works analogously in de Sitter 
space with our definition of the perturbative S-matrix. 

It thus remains only to address the more technical concern (i) . Dealing with the associated 
potential IR divergences (in this case, those not associated with self-energy corrections) is in 
fact the main focus of our work below. But the basic idea is simply that while most textbooks 
do use in/out perturbation theory to construct the Minkowski S-matrix, this choice is far 
from unique. The LSZ formalism allows one to extract the S-matrix from vacuum correlators, 
no matter how they have been computed [22]. In particular, one may define the Minkowski 
S-matrix by applying LSZ to time-ordered correlators computed using in/in perturbation 
theory via the closed-time-path formalism [23, 24]; see e.g. [25, 26] for reviews. To do so, 
one considers a path integral that begins at past infinity, runs up to future infinity, and then 
back to past infinity. One then inserts operators near past infinity (say, at the beginning 
of the closed-time path) and also at future infinity (in the middle of the closed-time path). 
Applying the LSZ formalism yields an S-matrix. The result is manifestly equal to the usual S- 
matrix, except that the out-vacuum has been replaced by the in-vacuum. But the positivity 
of the Minkowski conserved energy means that the in- and out- vacua agree, so that this 
closed-time-path definition does indeed reproduce the usual in/out S-matrix. Since in/in 
perturbation theory is well-defined in dS and leads to good asymptotic behavior (e.g., an 
analogue of cluster decomposition) [9, 8, 27, 10], the de Sitter analogue can lead to a good 
S-matrix. 

One complication is that in/in perturbation theory in de Sitter space is IR finite only 
when the closed time path begins and ends at a finite time. In particular, it is best behaved 
when the path begins and ends on a cosmological horizon so that it constructs correlators in 
the interacting Hartle-Hawking state (called \Q) below); see [28] for a proof valid for scalar 
theories with M 2 > 0. Thus the natural de Sitter analogue of the flat-space construction 
just outlined would be an S-matrix for the part of dS to the future of a cosmological horizon. 
This region is called the Poincare patch of dS by some authors and the expanding k = 
cosmological patch by others. While we believe that this object does in fact exist, we 
save its construction and analysis for future work. Instead, we construct an S-matrix for 
global de Sitter space below - roughly speaking by gluing together two copies of the Poincare 
patch construction just described, one for the expanding patch and one for the corresponding 
contracting patch to the past of the horizon - along the de Sitter horizon. More precisely, we 

2 That is, without resummation of the 1PI self-energy corrections to the propagator. 
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use an out/out-in/in closed time path which begins on a de Sitter horizon, travels backward 
to past infinity, returns to the de Sitter horizon and keeps on going to future infinity, and 
then finally returns to the de Sitter horizon once again; see figure 2. Operators are then 
inserted near past and future infinity and a de Sitter version of LSZ is used to extract an 
S-matrix. 



time 



Figure 2: The time contour used to construct our global de Sitter S-matrix. The surface S 
represents a cosmological horizon. See section 5 for details. 



This final LSZ-like step requires some care, as one must sort out further IR issues to 
find an appropriate dS analogue of the flat-space LSZ formalism. Some of these issues are 
already documented in the literature (see e.g. [29, 30, 31, 32]). The point is again that several 
procedures which are equivalent in flat space do not agree in dS, and that some seemingly 
natural dS generalizations lead to IR divergences for light (so-called complementary series) 
fields with masses < M 2 £ 2 < (D — l) 2 /4; here £ is the de Sitter length scale and D is 
the spacetime dimension. Nevertheless, we show in section 3.3 that a particular definition 
(based on extracting an S-matrix from poles in correlation functions) is free of unphysical 
IR divergences for generic 3 M 2 > 0. An alternate more explicit (but ultimately equivalent) 
procedure is described in appendix B. While we find it enlightening to understand this 
explicit structure, it quickly becomes cumbersome at higher orders in perturbation theory. 

Before proceeding, we should warn the reader of two further technical issues. The first 
is that our choice of time contour and the use of the associated Schwinger-Keldysh-like 
formalism means that perturbation theory involves several distinct types of vertices and 
propagators. This makes explicit computations more complicated than in flat space. It 
may well turn out that other computational techniques, such as analytic continuation from 
the Euclidean sphere or Euclidean AdS, will prove more efficient for computing scattering 
amplitudes. 

The second technical issue is that the lack of a positive-definite energy-like conserved 
quantity makes representation theory of the de Sitter group somewhat less powerful than its 
flat-space analogue. In particular, de Sitter representation theory alone does not guarantee 
the orthogonality of what we call distinct multi-particle asymptotic states. While such states 
are orthogonal for free theories, they cease to be orthogonal for general interacting theories. 
Even states with different particle numbers develop non-zero inner products. This, however, 
is not an obstacle to defining a useful S-matrix. We view it as essentially an accounting 
issue, though admittedly one that provides a sense in which our asymptotic states are not 

3 Our formulation fails for a measure zero subset of masses where self-energy corrections to the locations 
of poles are always large, even at small coupling g; see section 2.3. 
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'free.' A formal construction of orthonormal asymptotic states (e.g., by the Gram-Schmidt 
procedure) is always possible. However, such a construction obscures the physics, and so we 
will not emphasize it below. 

The remainder of the paper is organized as follows. We begin by reviewing basic proper- 
ties of de Sitter space and de Sitter QFTs in §2. We discuss the construction of asymptotic 
states in §3 and show that the resulting S-matrix has the desired properties in §4. The 
tools necessary to actually calculate S-matrix elements in Lorentz signature are provided in 
§5. Using these tools we compute the S-matrix of a model cubic theory of heavy fields to 
second order in §6. We then revisit our model theory in §7 allowing for fields of arbitrary 
positive mass. One interesting result is that, in contrast to the claims of [33, 34, 35], we 
find that complementary series particles generically decay. In all cases, we explicitly verify 
unitarity by showing that our S-matrix satisfies the de Sitter version of the optical theorem. 
We conclude with a summary and discussion of open issues in §8. 



2 Preliminaries 

We begin by briefly reviewing de Sitter space and some relevant aspects of de Sitter quantum 
field theory. 

2.1 de Sitter geometry 

The D-dimensional de Sitter manifold dSo may be defined as the single-sheet hyperboloid 
in an ambient (D + l)-dimensional Minkowski space: 

dS D = {X G M A1 \X-X = £ 2 }. (2.1) 

The line element of de Sitter space may written in the following convenient coordinates: 

7] G R. (2.2) 

Here £ is the de Sitter radius and denotes the line element on unit The 

time coordinate rj is related to the more familiar global de Sitter time coordinate t with 
g tt = — 1 via r] = sinh(t/^). In these coordinates the volume element is y/— g(x)d D x = 
£ D (1 + r] 2 Y D ~ 2 ^ 2 dr] d£Lu-\{x). The metric (2.2) describes spatial sections which are spheres 
S D ~ X of radius £ 2 {1 + rj 2 ). The conformal boundary of the chart consists of two disjoint 
spheres at rj —} ±oc denoted J? ± . 

At times it will be useful to use a second coordinate chart with line element 

r G K, (2.3) 

which also covers the entire manifold. Here dx 2 denotes the line element on M^ -1 . This chart 
has a coordinate singularity at r = corresponding to a preferred cosmological horizon T~L. In 



ds^ 



i + n 



;drj 2 + (l + rj 2 )dn 2 D _ 



ds 2 



= T 



-dr + dx 



zt2 
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the regions r < (r > 0) one may use instead the time coordinate A = ±r _1 to recover the 
more familiar Poincare coordinate chart (a.k.a. the expanding cosmological chart). Penrose 
diagrams which depict the above charts are shown at center (r-chart) and right (77-chart) of 
Fig. 1. 

While dS is maximally symmetric, it is important to remember that it has no globally 
timelike Killing vector field. In particular, neither d v nor d T are Killing vectors. As a result, 
global de Sitter does not possess a conserved quantity associated with flow in only timelike 
directions; in this sense there is no conserved "energy." For further details of de Sitter 
spacetime we refer the reader to [36, 37, 38]. 



2.2 Scalar fields on dS 

Free scalar QFT on a fixed de Sitter background is a well-understood subject with many 
good references (e.g., [37, 39, 27]). Below we focus only on establishing notation and on 
certain group theoretical aspects which will be useful later. At the end of this section, we 
comment briefly on interacting fields. 

Scalar fields on de Sitter are associated with representations of the de Sitter isometry 
group SO{D ) 1). For instance, the one-particle states of a free scalar field (j) a (x) described 
by the canonical Lagrangian 

1 M 2 
A>K] = --V^W.Or) - —<f>l(x), (2.4) 

form a unitary irreducible representation (UIR) of SO{D ) 1). We define the weight a via 

M 2 (a)f = -a(a + D-1). (2.5) 
The right-hand side is invariant under a —> — (a + D — 1); we choose a to be given by 



(D-l 



\2 



M 2 £ 2 



1/2 



(2.6) 



The UIRs of the de Sitter group may be classified as follows [40] : 

1. principal series: 

{ -^l<M 2 f, =► a = -^ + Vl pel, p>0, (2.7) 

2. complementary series: 

< mH 2 < {D ~ l) \ o e (-^^, 0) , (2.8) 



3. for some D, there is also a discrete series: 

M 2 f = -n{n + D - 1) for n € N , =>- a = n. (2.9) 
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"Heavy" fields belong to the principal series while extremely light fields with masses of order 
£~ 2 belong to the complementary series. These will suffice for our purposes as we limit our 
analysis below to fields with mass M 2 > 0. In contrast, UIRs of the discrete series are 
associated with M 2 = and certain tachyonic masses. It is useful to visualize the de Sitter 
representations in the complex a plane - see Fig. 3. 



Figure 3: Scalar de Sitter UIRs are depicted by the solid gray line in the complex a plane. For 
complementary series representations the weight a takes values along the negative real axis 
a G (— (D — 1)/2, 0) while for the principal series a takes complex values a = — (D — l)/2+zp, 
p > 0. For each series the 'conjugate weights' — (a + D — 1) are depicted with a dashed gray 
line. We denote by T p the Re a = — (D — l)/2 contour. Representations with a values to 
the left of T p are reducible; they may be represented as an integral over the principal series 
UIRs. Representations with Re a > — (D — l)/2 and Ima^O are not unitary. 

For a given mass M 2 > solutions to the Klein-Gordon (KG) equation (□— M 2 )0 cr (x) = 
are endowed with a positive-definite "Klein-Gordon norm." Working in the global chart, 
there is a natural basis of Klein-Gordon modes u a ^(x) which are distinguished by their 
angular momenta L on the spatial and are orthonormal with respect to the Klein- 

Gordon norm 



a 



i / (H7{x) \u a Zi (x)^ y ul l2 (x) 





(2.10) 



The Klein-Gordon modes may be written explicitly as 



(2.11) 



where Yg(x) are spherical harmonics on S 



parametrized by angular momenta L with 
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total angular momentum L and the functions f a L(v) ma y be written, e.g., as 



U(v) = N aL (l + V 2 ) L /\F l 
2 _ L _ (D _i )/2 



L-a,L + a + D-l:L + 



D 1 — irj 
2~ ' 2 



AT, 



r(^ + f) 



[T(L - a)T(L + a + D - 1)] 



1/2 



(2.12) 
(2.13) 



where 2 ^i(a,6;c;z) is the Gauss hypergeometric function [41]. The behavior of the mode 
functions in the asymptotic regions |?7| — ^ oo is given by [10] 

UM = (l + V 2 ) L/2 [^(^r-^ + K.^.!)^^)-^- 1 )^] [1 + 0( V - 2 ]] , 

M»(L-<t), (2.14) 



with K aL a coefficient whose value is 

2<r+(D-3)/2 



ah 



7T 



r a + 



£>- 1 



T(L - a) 



[T(L + a + D- 1)J 



1/2 



(2.15) 



It is important to note, however, that the Klein-Gordon current does not generally lead to 
a useful inner product between functions with support on multiple de Sitter representations. 
Integrals of the form 



i / 



u 



(2.16) 



aLi (x)^ u F(x) 

for general F{x) are complex, depend upon 7y , and may diverge as 1 770 1 —> 00. This includes 
the case where F(x) = u* r (x) is a Klein-Gordon mode corresponding to a mass a 2 7^ o. 
Indeed, the asymptotic behavior (2.14) shows that the expression 



lim 

?70— 7>±00 



(J2L2 



[x) 



V=V0 



(2.17) 



generally diverges for cti 7^ ct 2 , though it is finite for o\ = a 2 and is oscillatory for o\ 7^ a 2 
both in the principal series (in which case it converges to zero as a distribution). 

Free massive scalar fields in de Sitter have a unique normalizable maximally symmetric 
Hadamard state |0) known variously as the Hartle-Hawking (HH) state, the Euclidean vac- 
uum, and Bunch-Davies state [39, 42]. The multiplicity of names can be attributed to the 
multiplicity of ways to construct the state. The first two names come from the fact that the 
state, defined by its set of correlation functions, may be constructed by analytic continuation 
from the Euclidean section S D . Alternatively, one may construct this state in either global 
de Sitter or the Poincare chart by imposing adiabatic or "Bunch-Davies" vacuum conditions 
on the cosmological horizon. Within the static chart this state has yet another interpretation 
as the thermal state at the de Sitter temperature (it is the unique thermal state regular on 
the cosmological horizons). 
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For a scalar field 4> a (x) we denote the associated Wightman and time-ordered 2-point 
functions by 



W a {x u x 2 ) := (0| M^i)M^) |0) , (2.18) 

G a (x u x 2 ) := {Q\T4> rT (x l )4>a{x2) |0) . (2.19) 

They may of course be written in terms of the Klein-Gordon modes (2.11) as 

W a (x u x 2 ) = XX^ Xi Kl( X2 )' ( 2 - 20 ) 

L 

G a (x u x 2 ) = W (J (x u x 2 )0(r] 1 - r] 2 ) + W a (x 2 , £i)% 2 - Vi)- ( 2 -21) 



In the latter expression 0(rf) is the Heaviside step function taking values in {0, 1}. 

For interacting theories, perturbative corrections to such correlators often involve prod- 
ucts of several free 2-point functions W a or G a . Such computations can be simplified by 
making use of so-called linearization formulae, which express these products as weighted 
integrals of or over some contour in the complex /x-plane; see e.g. [8]. We will make 
use below of the slightly more complicated linearization formula 

^(J\L\ (^^(J 2 L 2 ^ 

) = ^CGC(L 1} L 2 ;£) (2n + D-l)p 

K ^ 

=: ^ f(2p + D - l)p 12 (/i, K)u^(x), (2.22) 

K ^ 

for the Klein-Gordon mode functions themselves (which are just de Sitter harmonics). To 
explain this formula, first recall that spherical harmonics obey their own linearization for- 
mula 

Y ti {x)Y t2 {x) = ^CGC(L 1 ,L 2 ;K)y i? (f), (2.23) 

K 

where CGC (...;...) denote the generalized Clebsch-Gordon coefficients of SO{D) [43]. 
Equation (2.22) then follows by using (2.11) and the analogous formula 

fa lL Mfa 2 LM= I (2fi + D-l)p aia2LlL2 (^K)f^ K (ri). (2.24) 

for the time-dependent parts of the modes u a ^. The kernel p ai a 2 L 1 L 2 (l Jj i K) ma y be computed 
the methods of appendix A of [8]. Indeed, up to normalization the p ai a 2 (p) f° un d there is 
our p ai a 2 L 1 L 2 (l Jj ^) evaluated at L\ = L 2 = K = 0. A useful Mellin-Barnes representation 
(A. 5) of p ai a 2 L!L 2 (p^ K) is given in our own appendix A, though we will not need the details 
here. 

In (2.24) the symbol J . . . denotes a contour integral in the complex p plane with 
measure dp/(2m). The contour is traversed from p = —ioo to p = +ioo within the strip 
Re(<7i + a 2 ) < Re p < 0. The representation (A. 5) can be used to show that p ai a 2 L!L 2 (p^ K) 
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is a meromorphic function of p within this strip and that it decays sufficiently rapidly as 
| Im/i| oo for the above integral to converge absolutely. At generic arguments 01, 02 the 
function p ai a 2 L 1 L 2 (l Jj i K) i n this strip has only simple poles at 

p = G\ + a 2 — 2n, \i — G\ — a 2 — D + 1 — 2n, 

/x = -<7i + a 2 - D + 1 - 2rc, /x = -a x - a 2 - 2D + 2 - 2n, for n G N ; (2.25) 

higher order poles arise at arguments <7i, a 2 where two or more of the above poles coalesce. 
These poles are required in order for the left-hand side of (2.24) to have the same asymptotic 
behavior in the regime \rj\ ^> L\ — 01, \v\ 3> L 2 — a 2 as the right-hand side. 4 It also useful 
to know that for <Ti, a 2 corresponding to positive mass-squared P cricr2 £ 1 £ 2 (M, ^0 satisfies 

so in particular p ai a 2 L 1 L 2 (pi K) G R for /x in the complementary series and generally complex 
for /x in the principal series. See appendix A for further comments. 

2.3 Interacting fields 

Finally, we review some relevant features of interacting massive scalar de Sitter QFTs. When 
the potential is bounded below, such theories admit a maximally symmetric state which 
is calculable at the level of perturbation theory, and for which perturbative correlators decay 
near past and future infinity. As a result, general states in the Hilbert space have the 
property that local correlators approach those of in the distant past and future. We 
say that \Q) is an attractor state for local operators [9, 27]. Since correlation functions in 

may be defined by Wick rotation from S D we refer to as the Hartle-Hawking state. 
However, these correlators may also be computed using the Lorentz-signature Schwinger- 
Keldysh formalism (see e.g. [28]). We delay a detailed discussion of theses techniques until 
§5.1 when such technical details become necessary. 

The structure of our S-matrix will be intimately tied to the de Sitter Kallen-Lehmann 
representation of the 2-point function (and to corresponding generalizations for higher cor- 
relators). Using the results of [8, 9, 27], it was argued in [44] that the 2-point function of any 
(perhaps composite) scalar operator $(x) in a general interacting theory may be written 

($(xi)$(x 2 )> := (fi|$(xi)$(x 2 )|fi> (2.27) 
p{p)W^x l) x 2 ) ) (2.28) 

where p(/x) is the Kallen-Lehmann weight of $ and / . . - denotes a contour integral in the 
complex /x plane (corresponding to complex operator weight /x) with measure dp /(2m). The 
integration contour is traversed from — ioo to +ioo to the left of the imaginary axis and to 

4 This can be verified by inserting the asymptotic expansion (2.14) for f^xiv) (2.22), closing the 
[i integration contour appropriately and using the Cauchy integral formula to equate the contour integrals 
with the sum of residues due to the poles (2.25). 
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X X 



Figure 4: An example contour of integration for the Lehmann-Kallen 2-point function. The 
contour is traversed from — ioo to +ioo in the strip to the left of the imaginary axis and 
to the right of any singularities (x's for our example) in the Kallen-Lehmann weight with 
Re a < 0. The particular example involves the theory discussed in section 6, which involves 
three species of scalar field. We have shown the poles in (<f>i(x)<f)i(y)) at 0{g 2 ) for a case 
where 0i, 02 lie in the principle series and 03 lies in the complementary series. 



the right of the singularities in p(/x) - see Fig. 4. See also [45, 46, 47, 34] for earlier related 
results. 

For example, in a free theory (2.28) is simply 



(2/i + L>- 1) 



W fl (x u x 2 ) = W a (x u x 2 ). 



(2.29) 



By deforming the contour of integration to lie along the UIRs one may convert the contour 
integral in (2.28) into an integral over positive mass-squared M 2 > (see Appendix A. 2 of 
[44]) but the form (2.28) is more useful for our purposes. Other 2-point functions (time- 
ordered, (anti-)symmetric, retarded/advanced) may be obtained by exchanging W M (xi,x 2 ) 
in the integrand of (2.29) with the appropriate free 2-point function. 

The singularities in the Kallen-Lehmann representation are particularly important. As 
emphasized in [8, 9, 27] they determine the large-?? behavior of the two-point function just 
as the behavior of flat-space correlators at large spacelike separations is governed by their 
singularities in the complex m 2 . Simple poles at UIRs contribute power law terms that 
appear in the corresponding free fields while higher poles contribute additional logarithms. 
More complicated singularities do not arise; in particular, there are no branch cuts 5 . This 
fact was shown in [9, 27], which analyzed the Mellin transforms of correlation functions in 
such theories in the perturbative expansion 6 . 

In perturbation theory about a free theory with weights c^, the following additional 
properties of p(/x) also follow from the Mellin-Barnes analysis of [9, 27]. These results will 



5 The familiar branch cuts of flat-space field theory arise in the limit i — >• 00 through the coalescence of 
an infinite family of poles. 

6 The Mellin transform of the 2-point function is closely related to the Kallen-Lehmann weight - see 
Appendix A.2 of [44]. 
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play key roles in constructing our S-matrix, showing that it is finite, and demonstrating the 
desired properties. 

i) Mass gap: All operators are "massive," by which we mean that their Kallen-Lehmann 
weight p(/x) is analytic on a strip in the complex \i plane e < Re/i < for some e < 0. 
As a result, there always exists an allowed path for the contour of integration. 

ii) Causality: The 2-point function of real scalar operators is real at spacelike separations. 
This requires that 

[/>(/*)]• = p(/i*). (2.30) 
Any singularity in the lower half-plane has an image in the upper half-plane. 

iii) Representation theory: The singularities of p{ji) are contained in the region 

Ru = {Re/x < - ^"^ U fie (-(£> - 1)/2,0)| . (2.31) 

This is the largest region consistent with scalar fields forming unitary representations 
of the de Sitter group. 

iv) Lattice of poles: For any scalar operator all poles in lie on the lattice built 
by taking linear combinations (with non-negative coefficients) of oi — 2n, —[(D — 1) + 
Gi] — 2n for non-negative integers n (n G No). Here oi are the weights oi of the free 
fields about which we perturb. 

We will call a theory exceptional if some weight c^, or some 'conjugate weight' &i = 
— [(D — 1) + Gj\ can be expressed as a non-trivial combination of the above generators 
with non-negative integer coefficients. By 'non-trivial,' we mean that the sum of the 
coefficients should be greater than 1; i.e., that they should not just be (1,0,0, . . . ) or 
some permutation thereof. Our definition of the de Sitter S-matrix will fail precisely for 
this (measure zero) set of exceptional theories. We will call a theory generic if it is not 
exceptional in the above sense. 

Note that for generic theories the bare weights Oi may be recovered given only the lattice 
of poles; i.e., given only the correlators computed to some perturbative order. 

v) Smooth self-energy: For generic theories, the one-particle irreducible (1PI) diagrams 
for two-point functions composite operators have no poles at the bare weights Oi or the 
bare conjugate weights c^. This may be shown by the same methods used in [9] to 
show that all diagrams (even if one-particle reducible) decay at large arguments at least 
as the two-point function of the lightest field involved. As a result, the corresponding 
self-energy II(/x) (see below) is smooth at \± — (Ji^Oi. 

Finally, when $(x) = (\) G {x) is a massive operator with bare mass M 2 (a) it can be 
convenient to write the Kallen-Lehmann form in terms of the self-energy II (/i) defined as 
usual by the sum of 1 PI 1 — >> 1 graphs: 
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In a strict perturbative expansion of the 2-point function about small couplings the denomi- 
nator of (2.32) must be expanded. Suppose the coupling constant is g\ then the perturbative 
expansion of II (/x) is 

oo 

nw^nw^), (2.33) 

71=1 

where 11^ (/i) is 0{g n \ and the expansion of the 2-point function is 



n<%) 



+ a + 1) 

n( 2 )(/i) / n«(/x) 



+ a + 1) \(fi-a)(fi + a + D- 1) 



. (2.34) 



Clearly then, order by order, each term in (2.33) inherits the five characteristics of p(/x) 
mentioned above. 

One may also consider the 2-point function obtained by truncating the sum (2.33) to some 
finite order TV and inserting the result into (2.32). While this object includes contributions 
from an infinite number of graphs at infinite orders in the strict perturbative expansion, its 
study is well-motivated when, say, the physics of interest is governed by the nature of the 
singularities (poles) in the Kallen-Lehmann weight. As in flat space, the self-energy can be 
used to compute perturbative shifts of these poles at each order. In generic theories, the fact 
that II(/i) is smooth at both the bare weights and the bare conjugate weights (see point (v) 
above) implies that these shifts are finite at each order. 

Correlation functions of three or more fields are then computed order-by-order via a 
diagrammatic expansion with the same general structure as various flat space constructions; 
see section 5.1 below. In particular, all higher-point correlation functions have a perturbative 
expansion analogous to (2.34) [27, 9]. Thus, for generic theories, the order of any pole at a i) 5^ 
in any n-point correlation function is controlled by the perturbative expansion in precisely 
the same manner as in flat space. 



3 Asymptotic states 

The quantum states involved in traditional scattering experiments may be described in 
terms of multiple widely-separated, non-interacting particles in the far past/future. While 
the construction of such particle states is trivial in free theories, we do not expect to be able 
to construct such states in interesting interacting field theories at finite separations. The best 
we can do is construct asymptotic particle states which enjoy a particle-state interpretation 
in the limits r\ —> ±oc. In the Heisenberg picture, inner products between past- and future- 
asymptotic states define the S-matrix (a more precise definition will be provided in §4). 

The construction of asymptotic particle states will present the main technical challenge 
below. We therefore devote this section to addressing this issue in detail. Our definition of 
asymptotic particle states is a straightforward generalization of definitions that appear in 
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standard QFT textbooks [48, 22, 49]. However, we warn the reader that one must use caution 
in applying the term 'non-interacting' to our asymptotic particle states. While our states 
transform in a simple way under the de Sitter group (and thus may be said to evolve freely 
in time) , for interacting theories the inner product on our space of asymptotic multi-particle 
states will differ from that associated with any free theory. 

3.1 General formalism 

The starting point for our construction is a (bosonic) local quantum field $(x) and the 
corresponding de Sitter Kallen-Lehmann spectral representation reviewed in section 2.3. As 
in flat space, we take the masses of the asymptotic particle states to be associated with 
values of a at which the integrand is singular. In particular, for interacting fields we take 
them to be given by the bare weights Oi defined by the lattice of poles; see point iv in section 
2.3. Precisely free fields have a Kallen-Lehmann representation supported only on a single 
dS UIR which defines the associated particles. 

Our goal is to use <&(x) to construct asymptotic particle states labeled by the quantum 
numbers of a de Sitter representation n = (a, L). We will give particular prescriptions for 
doing so below. For now, let us simply list the desired properties of such initial and final 
asymptotic states 

1. The states are normalizable: i/f{^\^)i/f < oo. 

2. The Hilbert-space is spanned by multi-particle states , i.e. any initial (final) state 
is a sum of states labeled by strings of quantum 7 numbers. So the states 1^)^ = 
|ni, n 2 , . . . , nk)i/f form a basis. 

3. Such basis states are symmetric with respect to their particle labels, i.e., |ni, n 2 , . . = 
1^2, ni,.. •)*//■ 

4. In the asymptotic past (future) the initial (final) asymptotic states transform as sym- 
metric tensor products under the action of the (diagonal) de Sitter group. This re- 
quirement is most easily understood in the Heisenberg picture. If U (g) is a de Sitter 
transformation with parameter g then a fc-particle asymptotic state transforms as 

U(g) \n u n 2 ,...,n k ) i/f = \gn u gn 2 , . . . , gn k ) i/f , gn { = (a, gLi). (3.1) 

The notation gL denotes an appropriate linear combination of angular momenta. The 
weight a is unchanged by de Sitter transformations because it denotes the eigenvalue 
of the quadratic Casimir of SO{D ) 1) which commutes with U(g). We say that states 
satisfying (3.1) transform covariantly under the the de Sitter group. We also require 
that (3.1) be unitary. 

7 Depending on the details of the theory more quantum numbers may be necessary to fully specify the 
state, but n = (a, L) will be sufficient for our purposes. 
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Each set of multi-particle states should include a vacuum state \v)i/f From criteria 1 
and 4 we conclude that this state must be both normalizable and invariant under the action 
of the de Sitter group: U(g) = \ v )i/f For free fields these requirements limit the 

possible choices of \v)^f to the so-called Mottola- Allen vacua [42, 39]. If we require that the 
short-distance structure be Hadamard, the free Hartle-Hawking state becomes the unique 
allowed choice . It is thus natural to take both \v) i ^ to be the Hartle-Hawking state \Q) 
in the interacting case as well. The fact that |Q) is an attractor state for local operators 
[8, 9, 27, 10] means that it is also a natural choice from the perspective of local physics in 
the asymptotic regions 8 . 

In the free case, the full space of states is readily constructed by acting with ladder 
operators 



|ni,ni, . . .,n k ) i/f := 



a ni a n 2 



l n k 



|o> 



where 



U at\ X . 



dZ u (x) k(i)V,«* £ (a;) 



i/f 



?7=const 



?7=const 



(3.2) 



(3-3) 



and as usual we also have 



S LK- 

(3-4) 

From the commutations relations (3.4) we readily verify that these particle states are normal- 
izable and symmetric, and from (3.3) we see that the ladder operators transform covariantly 
under the de Sitter group: U(g)a^ = a \ g ^ where gL denotes the covariantly transformed 
angular momenta. Unitarity of U(g) follows from the manifest de Sitter symmetry. Thus 
the states (3.2) satisfy our criteria for asymptotic particle states in both past and future. As 
an added bonus, this set of states also forms an orthonormal basis. Since 1 0)^ = \0)j = |0), 
the S-matrix is trivial. 

We will describe extensions of this construction to interacting fields in sections 3.2 and 
3.3 below. For now, we mention that the initial (final) states define two multi-particle Hilbert 
spaces T-Li/f and that for perturbation theory in small couplings we expect both Hi/f to be 
isomorphic to the entire Hilbert space T~L. So each set of multi-particle asymptotic states 
should provide a complete basis for T~L. We simply take this as an assumption, though our 
explicit verification of the optical theorem at low orders of perturbation theory in sections 6 
and 7 provides supporting evidence. 

We noted above that the states (3.2) are orthonormal for free fields. However, we have 
imposed no requirement on the inner product f (a\b) f between two final states (or i (a\b) i 



8 The authors [50, 51, 52, 53] have considered asymptotic states constructed out of the non-Hadamard 
Mottola- Allen vacua. Though they typically have divergent stress-tensors and do not reduce to the Minkowski 
asymptotic particle states in the flat-space limit [54], it has been argued that these asymptotic states may 
nevertheless play an important role in a dS/CFT correspondence. 
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between two initial states) for more general theories. In particular, because requirement (4) 
involves only the diagonal action of the de Sitter group (with all ^'s in (3.1) being the same), 
it is much weaker than requiring our asymptotic multi-particle states to be tensor products 
of asymptotic single-particle states, and this in principle allows the inner products f(a>\b)j 
to depend on the interactions. In Minkowski space, the continuous spectrum of the energy- 
momentum operator means that the associated conservation laws enforce orthogonality of 
distinct asymptotic particle states even in the interacting theory. But de Sitter space does 
not enjoy such simple conservation laws: there is no positive-definite conserved energy, and 
conservation of the spatial momenta is governed by the addition of angular momenta (i.e., 
the Clebsch-Gordon coefficients of SO(D)) which is less stringent than the addition of linear 
momenta. So long as \L\ — L 2 | < L 3 < L\ + L 2 , dS group theory alone does not require the 
vanishing of overlaps like ^xn^n^) { . Instead, we will be led to compute such inner products 
order by order in perturbation theory using appropriate Feynman-like diagrams. 



3.2 Interacting theories and IR divergences 



We now turn to the construction of asymptotic states for (massive) interacting fields; i.e., 
we seek a generalization of the LSZ formalism to de Sitter space. Unfortunately, the gener- 
alization of LSZ to de Sitter space is not unique. We describe one generalization below that 
is both conceptually and computationally straightforward. It turns out to suffice for explicit 
computations of heavy (principal series) fields but becomes IR divergent in the presence of 
light (complementary series) fields. Closely related IR divergences were noted previously in 
[29, 30, 31, 32]. After identifying the relevant issues below, we present a modified procedure 
in section 3.3 that is free of unphysical divergences. The procedure of §3.3 is also required 
to give a field-redefinition-invariant definition of the S-matrix even for heavy fields. 

Our first attempt at an LSZ-like prescription involves defining (now time-dependent) 
ladder operators in analogy with those of the free theory (3.3): 



We hope to construct initial (final) states via 



v 

(3.5) 



|ni,n 2 , . . . ) n k ) i 



i/f 



lim 

77—^00 



b=i 



(3.6) 



As mentioned above, asymptotic particle states will not generically remain orthonormal in 
the presence of interactions. But so long as the inner products are finite this need not be seen 
as an obstruction. Indeed, they may be computed order by order in perturbation theory as we 
discuss in sections 6 and 7 below. If desired, one may then to choose to construct orthonormal 
bases via e.g. the Gram-Schmidt procedure. The resulting orthonormal asymptotic states 
then enjoy all of the properties discussed in §3.1 except that they do not have well-defined 
particle numbers with respect to the field operator (j) a (x). 
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However, it turns out the inner products of these states can diverge even at lowest order 
in a coupling (i.e., at tree level) when complementary series fields are involved. In such 
cases the Klein-Gordon inner products used to define the ladder operators (3.5) may fail to 
have well-defined limits as rj — >> ±oc. Nothing is wrong with the theory per se: the (j) a (x) 
correlation functions are well-defined and for non-coincident configurations enjoy power-law 
decay as the operators are taken to large times \rj\ 3> 1. The problem is simply that the Klein- 
Gordon inner products contain a competing power-law growth dYiU a £ ~ (1 + r] 2 ) D / 2 r] a . The 
single potentially-dangerous term cancels precisely for free fields, but not with interactions. 



3.3 The R a projector 



The IR divergences found above appear to result from technical rather than physical issues. 
In Minkowski space, treating the KG inner product in a natural distributional sense usefully 
extracts the residue of the operator at a simple pole and allows us to define asymptotic states. 
When correlators of (\) G can be written as integrals over only principle series UIRs, the same 
is true in dS. But this procedure fails in dS when correlators receive contributions from 
complementary series UIRs. As already mentioned in §2.2, in de Sitter space the KG norm 
fails to converge on generic functions in the asymptotic regions, even when these functions 
decay exponentially. So it seems that the KG norm is simply not the right tool for extracting 
the complementary-series residues and thus the desired particle content from the correlation 
functions. 

We will resolve this technical issue by modifying (3.5) so as to remove contributions from 
any support away from the designated value of a. This may be accomplished by inserting 
the operator 

b 2 



lim exp 

b— >oo 



(D x - M 2 (a)) 2 



(3.7) 



into the definition of the ladder operators: 



a^ £ (ri) := -i J dE"(x) u al (x) V~ v R a ^(x) 



X 



(3. 



Note that R a annihilates smooth functions of a as well as poles elsewhere in the complex 
plane, though for any finite value of b it leaves the residues at a unchanged. Since all 
derivatives of R a with respect to a vanish at the pole, the higher order coefficients in the 
Laurent expansion about a are unchanged as well. Thus the multi-particle states defined 
by acting on \Q) with (3.8) are normalizable when the pole at a is simple, and these states 
agree with those built from (3.5) when all fields lie in the principal series. 

When the pole at a is not simple, it is important to point out that acting on \Q) with 
(3.8) still generally yields non-normalizable states. However, the remaining divergences are 
only logarithmic in 77, as opposed to the power-law divergences that arose for light fields in 
section 3.2. Similar logarithmic divergences are encountered in Minkowski space, where they 
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are associated with perturbative self-energy corrections that shift the original (bare) simple 
poles to new locations in the complex plane. The general structure of the 1PI expansion 
(2.34) shows that for generic theories (in which the self-energy II(/i) is smooth at the poles) 
the same is true in dS; i.e., the series (2.34) can be summed to yield (2.32) which has only 
simple poles near the original bare poles cr, — [a + D — 1] . While we will not work with the 
resummed series below, it is clear that logarithmic divergences of this type are physically 
meaningful. We therefore take them to be acceptable in an order-by-order computation of the 
perturbative S-matrix that does not use resummed propagators. We then say that, despite 
these divergences, our S-matrix "exists" for generic theories. We will comment further on 
the exceptional theories (in which II(/x) is not smooth at the poles) in section 8. 

An operator $ for which the naive LSZ expressions (3.5) converge up to logarithms will be 
called a "good operator". As we have discussed, examples include free fields of all masses, 
as well as general local operators in theories containing only principal series fields. The 
operator $(x) := R a (j) a (x) is a good operator in any theory. Some might object that this 
$(x) is a rather formal object, as the application of R a involves taking an infinite number 
of derivatives. But it is useful to note that there is an entire class of operators that are in 
some sense equivalent to $(x). Let us write $(x) = $(x) when the operators differ only by 
terms which decay faster than 0(r]~^ a+D ~ 1 " } ) 1 i.e. only up to terms which do not contribute 
to amplitudes of asymptotic states. Clearly any such $(x) is also a good operator. At a 
fixed order in perturbation theory it is straightforward (if tedious) to construct operators 
$(x) equivalent, in this sense, to R a (/) a (x) whose explicit expressions involve only a finite 
number of terms and contain only a finite number of derivatives of (\) a . This procedure is 
illustrated in appendix B. 



4 The S-matrix 

We have assumed that both the initial and final spaces of asymptotic states form complete 
bases for our Hilbert space. For generic theories we then define the S-matrix to be the 
operator which induces the corresponding change of basis. To be explicit, given bases \A)j 
and \B) i for the initial and final spaces of asymptotic states, we define the corresponding 
S-matrix Sba such that 

i^ = E^i^)/- ( 41 ) 

A 

For orthonormal bases, this gives 

S BA = f {A\B) i . (4.2) 

The essential steps in computing the S-matrix are the construction of correlators in the 
Hartle-Hawking state \Q) and the application to those correlators of the LSZ-like operations 
described in sections 3.2 and 3.3. Perturbatively, n-point functions in \Q) can be computed 
either by analytic continuation from the Euclidean D-sphere S D or equivalently (see e.g. 
[28]) by using Schwinger-Keldysh perturbation theory along a time-contour that begins and 
ends on some cosmological horizon. In order to insert operators near both I + and 7~, we take 
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the Schwinger-Keldysh contour to be as shown in figure 2. In the r-coordinates of section 
2, this may be described as the contour which for fixed x l starts at r = 0, runs backward 
to r = — oo, then runs forward to r = +oo, and finally returns to r = 0. We will use 
this scheme for our explicit computations in section 6 and 7. The associated diagrammatic 
expansion will be reviewed in section 5 below. 

However, regardless of the calculational scheme employed, the S-matrix enjoys the fol- 
lowing properties: 

1. The vacuum-to- vacuum amplitude is unity: By construction the initial and final 
vacuum states are explicitly the same state By definition this state is normalized 
to unity. 

2. Covariance under the de Sitter group: This follows from the covariance of the 
asymptotic states. E.g., 

f {A\B) i = f (A\l \B) f = f {A\U-\g)U{g) \B) f = f (gA\gB) i . (4.3) 

3. Behavior under CPT: Under the action of the CPT operator the S-matrix trans- 
forms as 

OS = S^G. (4.4) 

We note that for real fields like the scalar operators we consider section 6 and 7 the 
CPT operation is equivalent to the antipodal map A, so property (4.4) is simply a 
consequence of (4.3). 

4. Invariance under perturbative field-redefinitions: For simplicity, consider a the- 
ory defined by a single fundamental scalar field (j) a and consider the field 

= <t> a {x) + 0{x) (4.5) 

for some local operator &[x). We assume that &[x) may expanded in powers of c/) a . 
Since taking &{x) linear in the fundamental fields would induce only a trivial field 
renormalization (which might be matrix- valued if the theory involves multiple fields), 
we take this expansion to have no linear terms. As noted in section 2.3, the poles in the 
Kallen-Lehmann representation of G{x) lie on the same lattice as those of (j) a . They 
and are thus associated with the same bare weight a, and the same bare conjugate 
weight or, and define asymptotic particles of the same mass as <j) a . Since section 2.3 also 
noted there that for generic theories the operator &[x) has no poles at a, a themselves. 
It follows that R a Q = R a (j) a and thus that $(x) and (j) a define the same S-matrix. 

5. Unitarity: The physical content of S-matrix unitarity is implicit in our assumption 
that the initial and final spaces of asymptotic states each form complete bases for the 
same Heisenberg-picture Hilbert space, with the same norm. When expressed in terms 
of orthonormal bases, this assumption implies that the matrix Sba is unitary in the 
standard sense that S^S — 1 — SS^ . The unitarity of such an S is often phrased 
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in terms of the transition matrix T where 5=1 + iT . Then the unitarity of S is 
equivalent to the optical theorem 

2ImT=T t T. (4.6) 

While this remains an assumption at present, sections 6 and 7 will provide evidence 
in its favor by explicitly verifying the optical theorem at low orders of perturbation 
theory. 

6. Flat-space limit: Consider the flat-space limit £ — >► oc (with all masses and couplings 
held fixed in physical units). As mentioned above, the correlators used to define our 
S-matrix may obtained via analytic continuation from Euclidean signature; i.e., from 
the path integral on S D . The latter clearly approaches the usual flat-space result for 
large S D . In particular, the generating poles defined by the dS correlators approach 
the 1-particle poles of the Minkowski correlators. We noted above that the recipe 
(3.8) relates our S-matrix to the residues of such poles. Since this is analogous to 
the construction of the Minkowski S-matrix from residues of 1-particle poles, in the 
flat-space limit they must agree up to possible normalizations. The agreement of such 
normalizations will be verified by explicitly calculating the S-matrix for a 3 theory in 
section 6. We will find there that the de Sitter corrections are exponentially small in 
M 2 ^ 2 , and are thus nonperturbative in \ jt. 

5 Lorentz-signature perturbation theory 

The above sections give a general definition of our de Sitter S-matrix. In order to demonstrate 
the properties (1-6) from section 4, sections 6 and 7 below will explicitly calculate the S- 
matrix for cubic theories through 2nd order. We choose to use Lorentz-signature Schwinger- 
Keldysh perturbation theory [23, 24] (see also [55, 56, 13, 57]). Since we are unaware of a 
presentation in the literature of this technique in the global de Sitter chart we give a brief 
introduction in section 5.1 below 9 . Section 5.2 then states the associated diagrammatic rules 
for computing our scattering amplitudes. For now, these rules are valid only for sufficiently 
heavy fields such that there are no IR divergences. 

5.1 Schwinger-Keldysh correlation functions 

The generating functional for time-ordered Hartle-Hawking correlators in global de Sitter 
may be constructed as follows. Consider for simplicity a theory of a single scalar field (f)(x). 
We note that 

/ V 

9 The Lorentz signature discussion of [58] involves a similar 3-legged contour but uses the de Sitter "neck" 
r] = as the initial and final surface. In order to successfully produce the Hartle-Hawking correlators, this 
approach requires an independent computation of interacting correlators at r] = 0. This latter step is not 
required in our approach. See [28] for further comments on the use of spacelike initial surfaces. 
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where in the second equality we have inserted two complete sets of states {|/)} and {\p)} 
defined at a time in the far future and past, respectively, of an initial surface So- We 
take S to be a cosmological horizon so that regions to "future" / "past" correspond to ex- 
panding/contracting Poincare charts. The state \Q) may be defined by imposing adiabatic 
boundary conditions on S (see e.g. [28]). The expression (5.1) may then be understood 
as a path integral with a three-legged time contour of integration: reading (5.1) from right 
to left, the time contour begins at the horizon, is traversed into the far past, then to the 
far future, then back to the horizon. In the coordinates of (2.3), S is simply the r = 
hypersurface, the "future" / "past" are the regions r > 0/r < 0, and the time contour is 
r = 0— L — +L — for any sufficiently large L; see Fig. 2. 

By inserting classical sources J/, Jo, J p into the three legs of the path integral we obtain 
the generating functional for (j)(x) correlation functions with respect to 

Z[J f ,J ,J p ] := J2T,WMf\P^<PP)j, (5-2) 

/ v 



= /[Z%][2ty ][2ty p ]exp 



i {S[<l> f ] + J <f> f (x)J f (x)e,(x) 
+ i ( S[(f> ] + / 4> (x)J (x) 



i (s[(j>p\ + J (j) p (x)J p (x)9 p (x) 



(5.3) 

Here 9f(x) ) 9 p {x) are step functions with unit support in the future/past of the horizon, e.g. 
0(t), 9{—r). Time-ordered correlation functions of cj){x) are generated via 

<Q| T<f>( Xl ) . . . 4>{x n ) \Q) = - , . . . r Z[J p , J , Jf] 



iSJ (xi) iSJ Q (x ri 



(5.4) 



Jf=J o =Jp=0, 4>;- 



To generate diagrammatic rules one may, in the standard fashion, re-write the interacting 
part of the action S\ nt [4>] in terms of functional derivatives, then perform the Gaussian path 
integral over the fields. 10 The result is 



Z[J f , J , J p ] = exp (-iS^ 



with 

Z [J f , J p , J ] ■■= exp 




iSJ f 

1 

2 
1 



+ i3 



int 



s 



i8J 



— iS\ 



int 



5 



iSJ n 



Z [J f: J , J p ], (5.5) 



J (x)G(x,x)J (x) — -Jf(x)G*(x,x)Jf(x)6f(x)6f(x) 



J p (x)G*(x,x)J p (x)9p(x)9p(x) + J (x)W(x,x)J p (x)9 p (x) 



Jf(x)9f(x)W(x,x)Jp(x)9p(x) + Jf(x)0f(x)W(x,x)J o (x) 



(5.6) 



10 We assume the Lagrangian has the decomposition Co + £i nt with Co the canonical free field Lagrangian 
(2.4). 
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For correlation functions with arguments contained in a single Poincare chart the generating 
function (5.5) reduces to the more familiar generating functions of "in-in" or "out-out" 
closed-time-path perturbation theory [55, 56, 59, 60]. Likewise, for arguments contained 
within a static chart (5.5) reduces to the generating function of thermal field theory [26]. 

The generating function (5.5) differs from generating functional of "in-out" perturbation 
theory in that it has two additional path integrals (or "time contours"). These may be 
interpreted as providing the necessary corrections to (Q\ and \Q) respectively. Their origin 
lies in the fact that we have defined (Q\ and \Q) at an intermediate time rather than at 
the conformal boundaries rj = ±oc. In terms of Feynman graphs, these additional contours 
lead to extra vertices (which we label by / and p) in addition to those (which label o, for 
ordinary) arising a in-out" perturbation theory; i.e., we label a vertex /, o,p to indicate the 
term in (5.5) from which it arises. 

The resulting graphs can contain legs which are Wightman functions rather than time- 
ordered Green's functions. These have the important role of restricting the region of space- 
time integration involved in the construction of any correlator to the union of i) the past light 
cones of operators in the upper Poincare chart, and ii) the future light cones of operators in 
the lower Poincare chart. As a result, the perturbative construction of correlation functions 
via (5.5) is manifestly causal. 

5.2 Scattering amplitudes 

The above prescription for computing vacuum correlators leads directly to diagrammatic 
rules for scattering amplitudes. We work in the asymptotic particle bases and use the ladder 
operators (3.5), if necessary evaluating the states (3.6) at a finite cutoff r\ = ±r] . This 
allows us to deal with potential logarithmic divergences. The rj — >■ ±oc limits give rise to 
our S-matrix so long as the fields used in (5.4) are good operators in the sense of section 3.3. 
In particular, the rules below may be used if all fields involved are sufficiently heavy or if, 
in the case of complementary series fields, we have performed a field redefinition equivalent 
in the sense of section 3.3 to (j) — >■ R a (j)\ see appendix B for an explicit example. 

With the above understanding, we may interchange the order of the KG inner products 
and these spacetime integrals so as to first evaluate the KG inner products between mode 
functions and the Green's functions connected to external points, then take the rj ±oc 
limit, and finally evaluate the rest of the Feynman diagram. This is precisely what one does 
in the traditional Minkowski setting. 

Let us focus on the amplitudes of types i (a|fe) i , f(a\b) ^ or f(a\b) i for which we may use 
the time-ordered field correlators. Such amplitudes involve the p, o, and /-type vertices 
defined above according to the following prescription. First draw every Feynman diagram 
that contributes to an amplitude, paying no regard to different varieties of vertices. Then 
for each diagram: 

1. For every vertex connected directly to a final bra, include only p and o varieties, 
and replace the Green's function connecting to a final bra with where y is the 

integration variable. 
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2. 



For every vertex connected directly to a final ket, include only the / variety, and replace 
the Green's function connecting to a final bra with u n (y). See also rule 5 below. 



3. For every vertex connected directly to an initial bra, include only the p variety, and 
replace the Green's function connecting to a final bra with u* n (y). 

4. For every vertex connected directly to an initial ket, include only the / and o varieties, 
and replace the Green's function connecting to a final bra with u n (y). See also rule 5 
below. 

5. The diagram vanishes if it contains a vertex connected directly to both a final bra and 
final ket (in the computation of some ^(a\b)^) or to both an initial bra and an initial 
ket (in the computation of some i (a|6) i ). This removes any conflicts between rules 1 
and 2 and between rules 3 and 4. 

6. For every vertex that is not connected to an external state, include all varieties o, /, 
p. 

7. Include multiplicative factors of (ig) for every o vertex and — (ig) for every / and p 
vertex in the diagram. 

8. Each f/p vertex includes a theta function 9f/ p (y) where y is the integration variable. 

9. The Green's functions for lines between vertices are as follows: 

x vertex y vertex Green's function 









G„{x,y) 





f 


W a (y,x) 





V 


W a (x,y) 


f 


f 


G*(x,y) 


f 


V 


W a (x,y) 


V 


V 


Gl(x,y) 



These rules allow us to write scattering amplitudes as multiple spacetime integrals whose 
integrands involve products of KG modes, Green's functions, and the theta functions 9f/ p (y). 
These integrals are not easily computed in closed form but they they do reveal key properties 
of the amplitudes. 

A useful fool for analyzing amplitudes is the antipodal map A which maps a point x 
in global de Sitter to its antipodal point Ax. For a fixed coordinate chart (2.3), if x is in 
the past Poincare chart then Ax is in the future Poincare chart. A acts upon our basic 
ingredients as 

A : Y z {x) » {-l) L Y z {x), A:f aL ( V )^f a£ (-ri). (5.7) 
Complex conjugation acts upon our ingredients as 

Yi{x) = (-l) m Y £ (x), f; L (r ] ) = U(-v)- (5-8) 
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Here m is the quantum number associated with the azimuthal direction in the standard 
spherical coordinates. 11 Combining these results we obtain the relations 



= (-l) L+m K E (x), (5.9) 

W a (Ax,Ay) = W a (y,x) = [W a {x, y)}* , (5.10) 

G a (Ax,Ay) = G a (x,y), (5.11) 

9 p/f (Ax) = 6 f/p (x), (5.12) 

V-g(Ax) = y/-g(x). (5.13) 
The final equality follows from the fact that A is a de Sitter isometry. 

6 Example with heavy fields 

We now provide an explicit example of the construction of section 4. Consider a model 
theory of three massive scalars with a cubic interaction described by the Lagrangian 

3 

C[$\ = J2 £ 0[(f>i] + AntM +£c.t.[A (6.1) 

i=l 

£mt[<3 = gfofofaix), (6.2) 



2=1 



{Z^-l\ , V7M ,^, (Z Mi -l)Mi 



+ 0(g 3 ). (6.3) 



We work in perturbation theory with g<Cl, and we will examine S-matrix elements through 
0(g 2 ). The dS Hartle-Hawking correlators of this theory were studied previously [8]. The 
field and mass renormalization counterterms may be chosen to be = 1 + 0(g 2 ) and 
Zm % = 1 + 0{g 2 )] no other counterterms are necessary until 0(g 3 ). 
In this section we restrict the bare weights to the regime 

Re(ai + cj 2 + a s ) < —(D - 1). (6.4) 

This condition is automatically satisfied by principal series fields (for which Re a = — (D — 
1) /2), but it can also be satisfied when some or all of the fields are sufficiently heavy members 
of the complementary series. To 0{g 2 ) in the above theory, the condition (6.4) will suffice 
to qualify the fields as "heavy." In particular, the naive LSZ prescription (3.5) gives finite 
results and the diagrammatic rules of section 5.2 may be applied directly. We will revisit 
this theory in §7 in order to lift the restriction (6.4). 

Below, we use subscripts 1, 2, 3 to indicate which of the 3 flavors of scalar field in (6.1) are 
described by a given object. E.g. Ui,u[ denote wavefunctions of 0i-particles, n 2 ,ri2 denote 
occupation numbers for 2 -particles, etc. More primes will be added as needed. 



11 Our conventions for spherical harmonics are those of [61]. All we really need for our analysis is that 
Y^(x) oc e l7nip where ip is the azimuthal coordinate. 
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6.1 0(g) amplitudes 




Figure 5: Feynman graphs contributing to connected corrections at 0(g) and 0(g 2 ) in 
<?03^20i (oc) theory. We use solid lines for 0i(x), dashed lines for 2 and curly lines 
for 03 (x). In the top row, from left to right: the corrections (03(x3)0 2 (x2)0i(xi))^\ 
(02(x2)0i(xi)02(x2)0i(xi))^\ and (0i(x)0i(x))^ 2 \ On the bottom row are the mass and 
field renormalization counterterms present in (0i(x)0i(x))^ 2 \ There are also graphs of the 
same topology but with different permutations of the fields. 

Contributions to the S-matrix at 0{g) arise from the 3-point function (03(x 3 )02(x 2 )0i(xi)) 
Here and below we use parenthesized superscripts (e.g., to denote the order in g to which 
a quantity is calculated. We begin with the amplitude 

f (n s n 2 \n 1 )[ 1) = o term = ig J u^u^y), (6.5) 

where we abbreviate u\u<i . . .u* k (y) = u ni (y)u n2 (y) . . . (y). In the first equality we note, 
following the rules of §5.2, that only the V vertex contributes. The final integral expression 
is simply the de Sitter analogue of the familiar Minkowski result. 
The integral may be decomposed as 

/+oo i* 
drj(l + r?r-^f; 3L j; 2L J aiL M YZY£Y Ei (x). (6.6) 
-co J X 

The integral over x provides an SO{D) Clebsch-Gordon coefficient (CGC). These CGCs are 
always real, and this particular CGS is non-vanishing when the angular momenta satisfy the 
triangle inequalities [43]: 

\U - Lj\ <L k <L x + Lj, V{z, j, k} = {1, 2, 3}. (6.7) 

The remaining integral over r] is also real, as can be seen by changing integration variables 
77 — >■ — r] and using (5. 8). 12 Thus, the entire amplitude (6.5) is imaginary. Note that this 
integral converges only for masses satisfying (6.4). 

12 More generally one can use time reversal and parity to show that integrals of the form 

r k 1 / k 1 \ 

I~ / H^n^iy), Re $>* + $>;,■ <-(D-l), (6.8) 
J y i= i j= i \i=i j=i j 
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To obtain an explicit expression for (6.5) we will compute the integral 

h := / u 3 u 2 u{(y), (6.9) 



where we have taken a complex conjugate relative to (6.5) so that we may more directly 
employ the linearization formula (2.24) to replace the product u 3 u 2 with a weighted integral 
over a single mode function. 

We will also make use of two further tricks. The first involves the cutoff integral 

r r+vo r 

/ Ul u* 2 (y):= d V (l + V 2 Y D - 2 V 2 /«i«5(y). (6.10) 
Because m 1j2 satisfy the Klein-Gordon equation with mass M\ 2 respectively we have 

J u lU * 2 (y) = — 3 l _ M ^ J dE"(y) [u^uuKy)]^ - [m V^y)]^-^ . (6.11) 

This formula is valid for all complex values of o\ and o 2 apart from ^ = Lj + n or ^ = 
—Li — (D — 1) — n where n is a non-negative integer. 

The other trick is to note that when \q\ > 1 the functions f a L(v) ma y be split into two 
functions each of which solves the Klein-Gordon equation individually but contains only one 
asymptotic behavior. 13 We label the two behaviors "fast" and "slow": 

f*L{\v\>l) = KM + f* L (v), r*M = 0(rn, fl L (v) = 0( V -^ +D -V), (6 .i2) 

The explicit form of f^{rj) may be found using identities of the Gauss hypergeometric 
function, but we will only need the leading behaviors given in (2.14). Complex conjugation 
affects f s Jl{rj) differently depending on the series to which a belongs: 

\r*M\* = fld-ril [flMX = fl L (~v), a ec. series, 

lf* L (v)T = fL(-v), [flMY = f* L (-v), o € p. series. (6.14) 

Returning to (6.9), introducing the time cutoff \rj\ < ry , and using the linearization 
formula for u 2 u 3 (y) yields: 



h(H) = Y, [{2H + D- l)p 23 (/i, K) [ 



u uK u l(y) 



[ (2/i + D — l)p 23 (/i, L 1 ) ( u^uKy), (6.15) 



are real and vanish when either L tot = Yli=i Li or m tot = Yli=i m i 1S °dd- 

13 The function f a L(v) ma y he decomposed this way for |^| > 1, but the asymptotic behavior does not 
dominate the behavior until \rj\ ^> L — a [10]. 
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where the final step is just orthogonality of distinct spherical harmonics. The integration 
contour may be traversed anywhere in the strip Re (02 + 03) < /i < 0; for concreteness let 
the contour traverse Re \i = — e for infinitesimal positive e. We then integrate over y using 
(6.11) to obtain 



<*>--/{ 



(ji-<n)(ji + <n + D - l) 

<ffi"(y) [[^ £l V,<(y)] 



v=vo 



\ u »l^ v u i(y)\=-^ \- ( 6 - 16 ) 



Next we split f^iv) m t° f as t and slow parts as in (6.12). The f^* Ll {i]) contributes terms 
0(H ai+D - 1 - e ) while the f* L Jrj) contributes terms 0(H ai - £ ). Thus the latter, "fast-decay" 
terms vanish in the limit rj —¥ oo and we need only keep only the contributions of f^* L (rj). 
We now deform the integration contour to the left from Re/i = — e to Re/i = Re(a 2 + c 3 ) + e. 
From (6.4) we know 



Re(cr 2 + a 3 ) < Re -(<j x + D - 1) < Re a x < 0, 



(6.17) 



so the only poles we encounter in moving the contour to Re/i = Re(cr 2 + a 3 ) + e are the 
simple poles at /i — o\ and ji = —(ai + D — 1) (examine the denominator (6.16)). The 
remaining contour integral behaves like (9(^ e ( cr2+cr3 + cri )+ D ~ 1+e ) anc [ so vanishes as r]o — > oo. 
Thus 

h(Vo) = -{p2 3 K£i) J dE v (y) [[u s ^ v ul(y)] v=Vo - (*/)]„=_ 
+P2 3 (-(^i +D- l),Lt) J d^(y) ^u 1 ^ u ul{y)] v=ri0 - [u^ u ul(y)} v= 



-rjo 



(6.18) 



It is straight-forward to evaluate the remaining KG norms and take the limit rj — >> oo. We 
finally obtain 



I\ — — cot 



7T <7i + 



D — 1 



Pmivu L i) ~ P2s(-(<ri + D - 1), Li) 



(6.19) 



For u\ in the complementary series cot [n (a + ^^)] is real and so too is the second term in 
brackets, as P23(^ £ ^ ^1) £ ^- F° r a i i n the principal series cot [tt (a + ^y^)] is imaginary 
and so too is the term in brackets which is equal to 2i Im/^^i, So l\ is real in all 
cases. Buried within p2s(cri^L) is the SO{D) Clebsch-Gordon coefficient. 

The most salient feature of the result (6.19) is that the amplitude ^ (713712 = igh 
is generically nonzero whenever allowed by addition of angular momentum. This is true for 
any set of masses satisfying Re(oi + a 2 + 03) < (D — 1). This includes configurations in 
which some of the fields belong to the complementary series. The expression (6.19) contains 
P23(/^j L) evaluated in the region of the complex /x-plane for which ^23(^5 L) is analytic and 
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nontrivial; thus, it can vanish at most at isolated points. We emphasize this point in order to 
contrast with the statements of [33, 34] that complementary series particles can decay only 
in what we have called exceptional theories; this discrepancy will be explained in section 
8. Due to the detailed form of p ai a 2 L 1 L 2 (l Jj ^ ^0? the amplitude decays exponentially at large 
masses; see Appendix A. We have checked these results by direct numerical integration of 
(6.9). 

There is also a non- vanishing connected contribution to the amplitude 

(Q\n 3 n 2 n 1 )\ 1 ) = f + o terms = i / O p (y)u 3 u 2 u 1 (y) . (6.20) 

Jy 

This amplitude may be understood as both a transition amplitude for a 3 process and 
as a correction to the overlap between two initial states. As with the previous example (6.20) 
converges absolutely only for masses in the regime (6.4). The same physics may be accessed 
different ways: we may compute 

M^^V^-rJeM^M, (6.21) 

which is clearly the complex conjugate of (6.20), or we could compute the amplitudes 

f (n 3 n 2 n 1 \n) {1) = o + p term = ig / 6 f (y)u* 3 u* 2 ul(y), (6.22) 

(fi|ra 3 ra 2 ni)$i L) = / term = -ig J 9 f (y)u 3 u 2 u 1 (y), (6.23) 

which are also clearly complex conjugates. Moreover, by letting y —> Ay in any of these 
expressions we verify that 

.(n 3 n 2 ni|^) (1) = (^|n 3 n 2 ni)^ 1) , and f (n 3 n 2 n 1 \Q) W = {Vt\n 3 n 2 rii)f ) , (6.24) 

as is required by CPT. 

Evaluating these amplitudes explicitly is rather awkward because they are expressed as 
an integral of global KG modes over a single Poincare chart. We will be content to show that 
the amplitude is generically non-zero when L\ + L 2 + L 3 is even. Consider the imaginary 
part of the amplitude, 

2i Im <Q|n 3 n 2 ni>^ 1) = (Q\n s n 2 n 1 )[ 1) - (^n^n^p = ig J u 3 u 2 u r {y) =: igl 2 . (6.25) 

The integral I 2 may be evaluated by following the same procedure as was used to compute 
I\. The result is 



e-^cos [tt(^)] 



sin 



P 23(<TuL 1 )-p 23 (-(a 1 + D-l),L 1 ) . (6.26) 



Indeed, this expression is generically non-zero. We plot an example in Fig. 6 below. 

Note that all other corrections to either the initial or final inner products (e.g. j(ni|n 2 n 3 ) 
or j(ni|n 2 n 3 )^) vanish by rule 5 of section 5.2. 
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1 2 3 4 5 

M, 2 

Figure 6: The amplitude I 2 in D = 3 as a function of M\ with L\ = L\ = L 3 = 0, 
Mf = 2/^ 2 , and M| = 1.25/P. The continuous curve is plotted using (A. 5) while the dots 
represent direct numerical integration of (6.25). The amplitude is small, tends rapidly to 
zero as M\ increased, and approaches zero for M\ —> 0. Since M 2 and M 3 lie in the principal 
series, (6.25) converges for all M\ > - though the slow convergence near M\ « makes 
direct numerical evaluation of (6.25) difficult in that regime. 

6.2 0(g 2 ) state corrections 

We now compute corrections to the initial state inner products at 0(g 2 ). The corresponding 
final state corrections are then easily obtained by applying CPT. Since explicit expressions 
for amplitudes at 0(g 2 ) are both lengthy and time-consuming to attain, we will compute 
results only at the level needed to verify the optical theorem in section 6.3. Due to the way 
that the various terms enter the optical theorem, it will be convenient to delay computation 
of the actual S-matrix elements until section 6.3. 

At 0{g 2 ) there are two types of connected corrections to vacuum correlators functions: 
tree-level corrections to 4-point functions and 1-loop corrections to the 2-point functions 
(recall Fig. 5). Tree-level corrections to 4-point functions contribute to the connected part 
of the amplitudes 

i {n\n2n l 1 n' 2 \^l) c = pp term 

= (^#) 2 / / p (y 1 )e p (y 2 ) [uiU2(yi)u[u2(y 2 ) + u 1 u , 2 (y 1 )u , 1 u 2 (y 2 )} Gl(y U2 ), (6.27) 

as well as to amplitudes obtained by permuting the species. As we will see in section 
6.3, such vacuum-to-many amplitudes cancel out completely in any check of the optical 
theorem. There is also a tree-level correction to the connected part of the 2-particle amplitude 
i (^2 n il n 2^i)i c- I n principle both s- and u-channel graphs may contribute, though the u- 
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channel contributions vanish by rule 5 of section 5.2. Thus we compute only 

iK^|nin 2 )? ) _ channel =pf+po terms 

p (y)9 p (y)u^u / ^(y)u 1 u 2 (y)W 3 (y ) y) 



(ig) I °p(y) u i u 2 u a 3 K(y) 

y 



y J y 



K 



(6.28) 



These amplitudes are UV-finite because the SO{D) CGCs limit the range of the sum over 
K to 

min{|L' 2 - L[\ : \L 2 - L x \} <K< max{L' 2 + L[, L 2 + L x } . (6.29) 

The right-hand side can correct amplitudes both on and off the diagonal l^l' l 2 ' ® n ^ e 
diagonal the right-hand side is a sum of absolute values of integrals which were shown to be 
generally non-vanishing above. 

The 1-loop correction to the 2-point function (0i(x)0i(x))^ alters the normalization 
of the 1-particle states ^(ni|ni)^ , though conservation of angular momentum prevents it 
contributing to off-diagonal 1-particle amplitudes. The 1-loop correction is UV-divergent for 
D > 3, so we must include the standard perturbative counterterms that follow from (6.3). 
Only the field renormalization counterterm contributes to this amplitude: 



(2) 



= po + pf terms + (Z^ — 1) 

= -(wf 1 1 eMo P (yH(y)MyWs(y,y)w 2 (y,y) + (Z* - 1) 

Jy Jy 



EE 



L 3 



(ig) l_O p (y)u s u 2 u 1 (y) 
y 



(w) I O p (y)u 3 u 2 u 1 (y) 
y 



1). (6.30) 



This expression is IR-finite but the sum over momenta does not terminate and the expression 
contains an ultraviolet divergences for D > 6. We choose to cancel this divergence. All 
of this is very much analogous to the computation in Minkowski space. The 1-loop correction 

(0i(x)0i(x))^ 2 ^ also contributes a correction to the 2-to-vacuum amplitude ^1^0,^0^ , 

which will again cancel out entirely in any check of the optical theorem. 

There are of course analogous formulae for similar amplitudes obtained by permuting 
the flavors (1,2,3) above, as well as those between final states which may be obtained by 
using CPT symmetry and conjugation to relate final state amplitudes to the results above. 
For the remaining inner products between initial bras and initial kets which have yet to be 
mentioned, the 0(g 2 ) contributions vanish either manifestly (there are no relevant diagrams) 
or by rule 5 of section 5.2. In particular, we have 



>(2) 

he 



0, M^Kn^^O. 



(6.31) 
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6.3 The Optical Theorem 



We wish to explicitly verify the optical theorem, and thus the unitarity of our S-matrix, to 
order g 2 for our model theory. In terms of an orthonormal basis, the optical theorem takes 
the standard form -v_ 

(6.32) 



2Re f (B\A) i = y £ i (B\C) ff (C\A) i 



c 



where ^(B\A) i := ^(B\A) i — Sab- We may restrict attention to amplitudes involving at most 
4-particles as only these have fully-connected contributions. 

To verify (6.32), we first construct orthonormal bases from our multi-particle asymptotic 
states. At O(g) the 1- and 2-particle states require no modification; we need only remove 
the overlap of 3-particle states with the vacuum, as well as overlaps between higher-number 
particle states due to disconnected contributions. We therefore define 



nsn 2 n 1 )? N = \n 3 n 2 n 1 ) i - (^|n 3 n 2 ni)^ ij + 0(g 2 ). 



(i) 



(6.33) 



The only non-zero connected contributions to the S-matrix at 0(g) are the amplitudes 
/ ( n ?> n 2 \ n i)i which were computed in section 6.1. Recall that these amplitudes are imaginary, 
so that (6.32) is satisfied at 0(g). 

The vacuum-to-three amplitude in the orthonormal bases vanishes by construction: 



°J(n 3 n 2 n 1 \n)f N{1) 



0. 



(6.34) 



Quite generally, converting in this way from a particle basis to an orthonormal basis removes 
contributions from any diagram which is attached to only the final bra or initial ket (in the 
particle basis). It is natural to think of such diagrams as computing corrections to the 
relevant bra or ket, rather than representing genuine particle scattering processes, so it is 
not surprising that these diagrams do not ultimately contribute to the S-matrix despite being 



non-zero. 



At 0{g 2 ) ) transforming to orthonormal bases once again removes overlaps between parti- 
cle states and the vacuum; the transition also removes overlaps between non-identical pairs 
of 2-particle states and corrects normalizations. Below we use the orthonormal states 



I \ON I v 

\ n i)i = K) 



1 / \° N 1 f \ 

\n 1 n 1 ) i = l^ni) 



0(g 3 ), 



1 - ^K^iK^i 



,( 2 ) 



+ 0(g 3 ), 



\n 2 ni)° N = \n 2 n l ) l \ n 'i n 2)i M'M^df* + °{^)- 



(6.35) 
(6.36) 
(6.37) 



Returning to (6.32), we begin with the — >■ scattering process. Again, angular 
momentum conservation sets the off-diagonal terms to zero so that the non-trivial part is 
just 



°7< 



(nMf*® = f ( ni \ ni )\ 2) - I ifafa)® + finMf = (graph) + (c.t.) (6.38) 
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Here (graph) denotes the contribution from the 1-loop diagram while (c.t.) denotes contri- 
butions from the mass- and field-renormalization counterterms. Explicitly these are: 

(graph) = of + oo + po + pf — - [po + pf + of + pf] terms 

= oo H — [of + po] terms 

= (ig) 2 ! ! u* 1 (y)u 1 (y)M(y,y), (6.39) 

where 

M(y),y := G 3 (y,y)G 2 (y,y) - \w 3 (y ,y)W 2 (y ,y) (d f (y) + 6 p (y)) , (6-40) 

and 

(c.t.) = o term = i(Z Ml + -2) f u\ Ul {y). (6.41) 

Note that the integrand in (6.41) is positive-definite and behaves at large \q\ 3> 1 like 
r] 2ai so that the integral is IR-divergent for all values of o\. This divergence turns out to 
cancel a similar divergence that appears in (graph) under precisely the same conditions that 
UV counterterms are necessary. There is also a logarithmic IR divergence proportional to 
n^(cri) as expected from the second term in (2.34), which we therefore regard as physical. 
Recall that a similar divergence in Minkowski space is associated with a finite decay rate per 
unit time ("Fermi's golden rule") and plays a key role in satisfying the optical theorem. 
To verify our dS the optical theorem we compute 

— 2 Re ^(^iI^i)^ = —(graph) — (graph)* 

= -(ig) 2 [ f u\(y) Ul (y) [M(y, y) + M*(Ay, Ay)] . (6.42) 
Jy Jy 

In the second equality we have used (5.9). From (5.11) and (5.10) it follows that 

M(y,y) + M* (Ay, Ay) = G 3 (y,y)G 2 (y,y) + G* 3 (y,y)G* 2 (y,y) - W 3 (y,y)W 2 (y,y) 

= W 3 (y,y)W 2 (y,y). (6.43) 

Inserting (6.43) into (6.42) we obtain 

-2Re°> 1 |n 1 )f iV(2) = -(ig) 2 [ [ u* 1 (y)u 1 (y)W 3 (y,y)W 2 (y,y) 



y J y 



EE fa) / u ± u *2 u Uy) 

L2 L3 



EEl^k)^ .(6.44) 

L2 Lz 



The final equality is precisely what is required by the optical theorem. 
Next we examine the (f>2<fii — > 4>24>i amplitude 



v (2) 1 r 



J , (nir^|rein 2 )^ vw = /K^K?^}-^ - ^ /K^K^)/ + jK^l^i^)! • (6.45) 
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Figure 7: Corrections to 2 0i —> 2 0i scattering at 0(g 2 ). From left to right are the 
s-channel, the u-channel, and two disconnected 1-loop corrections. 



This expression naturally splits into a connected part f(n f 1 n , 2 \n\n2)fc which receives tree- 
level corrections and a disconnected part that receives 1-loop corrections; see Fig. 7. The 
connected part may further be split into s-channel and u-channel contributions with the 



former including the state corrections — | 
agrammatic rules it is straightforward to compute 



.(2) 



. Using our di- 



°^(n;n / 2 |n 1 n 2 )^I^ annel = +oo + of + po + pf - ^ [pf +po + pf + of] 

= +oo+ ^(of+po) 



\ON(2) 



with 



(ig) 2 / u , ^u , ;(y)u 1 u 2 (y)M{y ) y) ) 



M(y,y) := G 3 (y,y) - -W 3 (y,y) (6 p (y) + 6 f {y)) 



(6.46) 



(6.47) 



To verify the optical theorem we compute 

-2Re°J{n[n' 2 \n 1 n 2 )f^ aQnel = -(igf ! f u'*u'*(y) Ul u 2 (y) [M(y,y) + M*(Ay, 

Jy Jy 



(6.48) 

where in the second term we have once again used (5.9). From (5.11) and (5.10) it follows 
that 

M(y,y) +M*(Ay,y) = G 3 (y,y) + G* 3 (y,y) - W 3 (y,y) = W 3 (y,y), (6.49) 
and inserting this into (6.48) we obtain 

-2BB°J(n[n' 2 \n 1 n 2 )f^ = ~(^) 2 /_ / u^u , ;(y)u 1 u 2 (y)W 3 (y,y) 

Jy Jy 



L 3 



(ig) / u 3 u' 2 *u'^(y) 



(ig) / U3tt 2 «i(y) 



M'2 n 'i\ n ^)f f {nz\n2ni)f ] , 



(6.50) 
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where the summand vanishes unless 

max{|L 2 - L[\, \L 2 - L x \} < L 3 < min{L 2 + L[, L 2 + L{\ . (6.51) 

This is indeed a piece of the necessary result. To complete the desired check of unitarity, we 
must show that the u-channel amplitude combines with the 1-loop contributions to yield 

od OfM / I \OJV(2) 

-2 tie f {n 1 n 2 \n 1 n 2 ) i u _ channe i+ 1 _i oop 

= E E E MWiKnWD? fWn'Xlnin^ 

fill f ll j ll 
L, 1 L x 

+ E E E ^^iK^f ^nln'^n^ . (6.52) 

f III fll Til 
^2 ^3 ^2 

It is easy to repeat the above steps for the u-channel, as there is only one term. Doing so 
verifies that it contributes the part of the right-hand side of (6.52) arising from a cut tree 
diagram. Our earlier result for — 2 Re f(N\\ N\)f^ shows that it indeed provides the part 
of the right-hand side of corresponding to cut loop diagrams. Thus the optical theorem is 
satisfied. 



7 Light fields 

Section 6 restricted study of the model theory (6.1)-(6.41) to bare weights satisfying a\ + 
^2 + ^3 < _ (D — 1). We now extend the treatment of this model to arbitrarily small positive 
masses using the R a projector of section 3.3. As before, our goal is to verify the optical 
theorem to 0{g 2 ). A more explicit but ultimately equivalent method of computing the S- 
matrix for light fields is described in appendix B, which studies the same model theory to 

0(g 2 )- 

Let us begin by analyzing amplitudes at O(g). It is important to note that the actions 
of R a (x) and time-ordering do not commute, e.g.: 

Rt^R^x) (7>i(a;)0i(x)) = (T [R^^R^x)}) + (contact). (7.1) 

While these contact terms ultimately do not affect scattering amplitudes, they do complicate 
the computations. We find it simpler to first apply the Ri(x) operators on the Wightman 
functions of the theory, andjhen to use the diagrammatic rules of §5.2 to compute ampli- 
tudes. It is useful to define R\(x) through 

R a {x) = l + R a (x){D x -M 2 {a)), (7.2) 

so that 

Rz{x Z )R 2 {x2)Rl{x 1 ) (0 3 (£3)02(£2)0l(£l)) (1) 

= (4>z{x3)4>2{x2)4>i{xi)) {l) -gRz{x) (02^1 (^3)^2(^2)^1 (^i)) (0) 

-gR 2 (x) (0 3 (x 3 )030i (^2)0i(^i)) (O) - gRi(x) (0 3 (:c3)02(£2)0 3 02(£i)) (O) 

= (03(^3)02(^2)01 (^i)) (1) - gR 3 (x)W 2 (x 3 , X 2 )Wi(x 3 ,Xi) 

-gR 2 (x)W 3 (x 3 , x 2 )W 1 (x 2 , xx) - gRi(x)W 3 (x 3 , x 1 )W 2 (x 2 , xi), (7.3) 
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where we abbreviate Ri := R ai . 

Let us compute the amplitude j(n 3 n 2 |rii)^. We may time order the first term in (7.3) 
and use our diagrammatic rules. The contributions of the remaining terms in (7.3) are easily 
evaluated because they are free theory correlators: 



(n 3 n 2 \ni 



ig 



ulu* 2 ui(y) - g I dY? 



u 3 u 2 ) 



9 
2 



dZ u 



(i? 3 M2tXl)V V U* 3 



+ (2^3) 



-■no 



(7.4) 



Note that we have temporarily implemented a cutoff \rj\ < rj , and that we have explic- 
itly symmetrized across the final state operators for convenience. We may clean up this 
expression by writing the second term as an integral over the regulated space: 



(n 3 n 2 \n 1 



ig 



u* 3 u* 2Ul (y) + ig Ul (D- M^R^u® 



y,vo 



9 
'2 



/ dZ u \(R 3 u* 2 u 1 )V'„u* 3 \ + (2 o 3) +g \ dYT L ^ u {Riu* 3 u* 2 ) 
J L \+ m J J L 

ig j (R lU * 3 u* 2 ) Ul (y) - 9 - \ dY? V%u* 2 u x ) V„d + (2 o 3) 

J v ,rm Z J 1 J+7 ?° 



+9 / dZ v 



Ui 



^v(Riu* 3 u* 2 ) 



+V0 



We need to take a moment to analyze integrals of the form 



(RyuluDu^y). 



(7.5) 



(7.6) 



This integral is finite for all positive masses. When the weights satisfy Re(o"i + a 2 + cr 3 ) < 
— (D — 1) we may integrate by parts to make R\(y) act on u\(y)\ noting that R\{x)u\{x) = 
ui(x) it follows that 



h = / {RiuluDu^y) = / ulu^ur 

J V J V 



(y) = I u for Re(<7i + a 2 + a 3 ) < — (D - 1). (7.7) 



Using the linearization formula for utu^ we may compute 7 3 in precisely the same way we 
computed I\ in §6.1. The only difference is that now the R\{x) annihilates any contributions 
from poles in the complex \i plane other than the two at /i = a and fi — — (a + D — 1). Thus 
for all positive masses 7 3 is given by 



h = (RHS of eq. (6.19)) VM^ > 0. 



(7.8) 



It follows that in fact 

h 



f {RyuluDuM = I {R 2 u* 3 u l )u* 2 {y) = [ {R z u* 2 u\)u* 3 {y) . (7.9) 
y Jy Jy 
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All three expressions are analytic (up to poles at the exceptional theories). So (7.9) is clear 
from the fact that the integrals agree for Re(oi + 02 + 03) < — (D — 1) when the actions of 
i?2, R3 are trivial. 

We now return to the amplitude (7.5). Noting that 
(Riulu* 2 )ui(y) = / ulu* 2 ui{y) 

+ [ d^y L V^(^i^)l - / dYy LV.^i^^)] (7.10) 

y L J -770 7 L J +770 

and using (7.9), one sees that the surface integrals in (7.5) cancel each other completely so 
that the final expression for the amplitude is given by ig times (6.19) for all > 0. We see 
that the the amplitudes are analytic in the weights c^, though with poles at the exceptional 
theories 14 defined in point iv of section 2.3. As in section 6, the 0(g) S-matrix element 
f (^3^2 l^i)^ is imaginary in agreement with the optical theorem. 

A similar analysis (with corresponding results) can be performed for the remaining am- 
plitudes at 0(g). In particular, for all positive masses we have 

2ilm(n\n 3 n 2 n 1 ) i = (ig) (R.H.S. of (6.26)) V M- 2 > 0, (7.11) 

see Fig. 6 for an example. The construction of an orthonormal basis also proceeds just as 
for heavy fields and so again has no effect on the 1- and 2-particle states. 

We now proceed to 0(g 2 ). Once the R a operators are in place the calculation is similar 
to that described in section 6.3 for heavy fields. We will thus proceed rather quickly. 

The inner products which vanish for heavy fields continue to do so for lights fields as well, 
e -S- i( n 2 n i n 2|^i)ic = 0, so it is natural to again define orthonormal states via the analogues 
of (6.35). Since the quantities of interest for the optical theorem are the real parts of the 
0(g 2 ) S-matrix amplitudes f(A\B)f \ we focus only on these objects below. 

We begin with the 1—^1 S-matrix element °^{n\\n^° N ^ 2 \ The real part may be 
expressed as 

- 2Re°> 1 |n 1 )? iV(2) = MMi)? + f^i^f - - ^n^f . (7.12) 

Integrating by parts and using the fact that the (j)\(x) 2-point function is de Sitter-invariant, 
one may recast the right-hand side as 

-2Re°> 1 |n 1 )f iV(2 )= f u 1 (y)(D y — Mf) U{y){U ¥ - M^R^R^y) (faffiMv))® • 

Jy Jy 

(7.13) 

We may freely commute the Klein-Gordon operators through the R\ projectors and let 
them act on the correlation function. The KG operators reduce the correlator to an O(g ) 



14 These poles precisely match those found in [33, 34]. The coefficients in (B.7) also have poles at these 
theories. 
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correlator: 



-2Re°> 1 |n 1 )f iV(2) = +g 2 f Ul (y) [_u* (y) Ri(v) Ri(V) {^Uv^Uv))^ 

= +i (My) f «'®)Ri(v)RiW>[WMv)w 2 ($,y)] 



L2 L3 



id / {Riu* 3 u* 2 )ui(y) 



(7.14) 



The right-hand side is precisely what is required by the optical theorem. 

We conclude this section by verifying the optical theorem for the 2—^2 S-matrix element 
J(n , 2 n , 1 \n2n1 )f N • We focus on the s-channel contribution as this is the most non-trivial; 
the u-channel contribution to the optical theorem involves fewer diagrams, and the discon- 
nected contributions follow from (7.14). The real part of °^{n f 2 n f 1 \ri2ni)f N ^ may be written 

- 2Re°^(n / 2 n / 1 |n 2 ni)f 7V(2) = j [ [ [ [ Iu[%x 1 )u^(x2)u 2 (x2)u 1 (x 1 ) 

- Mt){U^ - Ml){U X2 - M 2 2 )(D X1 - Ml) 
R 2 (x 2 )Ri(xi)R 2 (xi)Ri(xi) 

{{<h&), Mxi)H<hW, 0i(^i)}) (2) )■ (7.15) 



We have explicitly symmetrized the states for convenience. Because the correlator is 0{g 2 ) ) 
only two Ri(x) operators act non-trivially on the correlator. Let us define the a s-channel 
part" to be the part for which 2 (x 2 ) and 0i(xi) are connected by a vertex. Then the 
s-channel contribution to (7.15) is 

-2Re ^(nXI^2n^^ [ [ [ { ^ / 1 *(xi)4*(x 2 )^ 2 (x 2 )^i(xi) 

(□* a - M 2 2 )(D 3?1 - Ml){U X2 - M 2 2 ){U Xl - Ml) 

[R 2 (x 2 )R 2 (x 2 ) + R 2 (x 2 )Ri(x 1 ) + R 1 (x 1 )R 2 (x 2 ) + R^x^R^xx)] 

({<fa(x 2 ), (f>l(xi)}{fo( X 2)i M X l)})s- channel 

(7.16) 



Once again we may freely commute the Klein-Gordon operators past the Ri(x) operators 
and let them act on the correlator. Note that 

(Q, 2 - Ml)(U^ - Ml)(U X2 - M 2 2 )(U X1 - Ml) <{0 2 (x 2 ), ^fa)}, {02(x 2 )0i(x 1 )})f_ ) channel 
= +9 2 S(x 2 ,x 1 )S{x 2 ,x 1 ) [W 3 (x 2 ,x 2 ) + W 3 (x 2 , xt) + W 3 (x u x 2 ) + W 3 (x 1 ,x 1 )} (7.17) 

Upon inserting (7.17) into (7.16) we obtain a number of terms. In each term we may use the 
delta functions to integrate over one of the Xi and one of the X{. For each term we expand 
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the Wight man function W${xi,Xj) in a sum over harmonics. Cleaning up a bit and using 
(7.9) yields the desired result: 



o -T) ON i i i I v 

— 2 tie f (n 2 n 1 |n2ni) 



OTV(2) 



# 2 (Rsu f 2 u[(y)) I {Rzu 2 u 1 {y))W z %y) 



£ 3 



(ifif) / (R3U 2 u 1 )u* 3 (y) 



(7.18) 



which explicitly verifies our optical theorem. 



8 Discussion 

We have introduced a perturbative S-matrix for interacting theories of massive (M 2 > 0) 
scalar fields in global de Sitter space. Our S-matrix reduces in the flat-space limit to the 
usual Minkowski S-matrix, is invariant under perturbative field redefinitions, and is uni- 
tary in the sense stated in section 4. In particular, we have verified the associated optical 
theorem through order g 2 in a gcj) 3 theory by explicit calculation. While there are certain 
qualitative differences from Minkowski space (lack of branch cuts, non-orthogonality of dis- 
tinct asymptotic multi-particle states) , we expect our S-matrix to be a useful tool for better 
understanding interacting de Sitter quantum fields. 

As expected from the lack of a positive-definite energy-like conserved quantity, we find 
that all particles generically decay for any M 2 > 0; i.e., all particles in dS are unstable. 
Let us pause briefly to explain the discrepancy with [33, 34, 35], which claim that general 
complementary series particles can decay only in theories that are exceptional in the sense 
of section 2.3. Although refs. [33, 34, 35] do not define a full S-matrix, they study what 
they call the inclusive decay rate V per unit time at 0{g 2 ). They compute the result in the 
principle series (where the IR issues are straightforward) and analytically continue to the 
complementary series. Their T is defined by generalizing certain formulae from Minkowski 
space, and in the principle series it agrees precisely with the right-hand-side of our 0(g 2 ) 
optical theorem (6.44). In the principle series it is thus —2 Re °^(ni\ni)f N ^ . The technical 

discrepancy with our work then arises because while the full amplitude °J(rii\rii)f N ^ is 
analytic in our formalism, the real part alone fails to be analytic outside the principle series. 
The issue is analogous to attempting to compute Ree^ = 1/e at z = i by analytically 
continuing from the real axis the function cos(x) = Kee lx and obtaining cos(i) = e + e _1 . 
The forthcoming work [62] will comment further on this issue. 

In Minkowski space, unstable particles are ultimately unsuitable as asymptotic states. 
In particular, this issue can be seen in perturbation theory once self-energy corrections are 
included in external lines representing the unstable particles. Before resumming 1PI graphs, 
the correlation functions develop higher order poles that lead to IR divergences in the S- 
matrix. After resumming 1PI graphs and interpreting the result as a self-energy correction, 
the relevant poles move off the real axis and the associated wavefunctions are no longer 
normalizable. 
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We find similar behavior for unstable particles in the principal series of de Sitter rep- 
resentations; i.e., with M 2 £ 2 > (D — l) 2 /4. Thus, a theory with only principal-series par- 
ticles apparently does not have a well-defined S-matrix outside of perturbation theory. In 
Minkowski space, the notion of a generalized S-matrix with unstable particles as additional 
asymptotic states has been sometimes been considered (see, e.g., section 4.9 of [63]). In 
de Sitter space with principle-series fields only, such a formalism could be used to define a 
generalized S-matrix despite the lack of stable particles. 

The behavior of complementary series fields is more subtle. As always, the optical the- 
orem relates the l-to-2 decay amplitude to the imaginary part of the loop-corrected 1-to-l 
amplitude. However, as will be explained in detail in [62], the fact that the wavefunctions 
do not oscillate near infinity means that the self-energy remains real; i.e., the imaginary part 
of the 1-to-l amplitude is not then given by the imaginary part of the self-energy, and the 
poles remain on the real axis. This raises the possibility that one might be able to define a 
finite S-matrix using normalizable asymptotic states built out of complementary-series fields, 
though this seems far from straightforward; see [62] for additional comments. 

With the above understanding, our S-matrix is free of unphysical IR divergences for 
generic masses. However, extra IR divergences do arise when the self-energy n is singular at 
some generating pole; i.e., for the exceptional theories described in section 2.3. As a result, 
the locations of the poles fail to admit a good perturbative expansion. In this sense, such 
exceptional theories are always strongly coupled even for small g. Our exceptional theories 
coincide precisely with the special cases previously identified in [33, 34, 35]. 

There are of course many possible generalizations of our construction above. For example, 
we worked here in global de Sitter space. We expect that one could use the future half of 
our Schwinger-Keldysh contour (see figure 2) to define a similar S-matrix in the expanding 
Poincare patch, though we have not explored this case in detail. In addition, we have thus 
far restricted attention to scalars with M 2 > 0. Massless minimally coupled (M 2 = 0) 
scalars do not have a good dS Hartle-Hawking state in perturbation theory (see e.g. [64]), 
but massless higher spin fields are better behaved. In particular, both the free Maxwell field 
and graviton admit de Sitter-invariant vacuum states with good analytic continuations to 
the Euclidean sphere (i.e., Hartle-Hawking states). [65, 66, 67] One may therefore hope to 
extend our S-matrix to such cases. So long as the associated correlation functions decay 
sufficiently rapidly 15 , our use of Klein-Gordon inner products (and the invariance of such 
inner products under linearized gauge transformations) suggests that the final result will 
be gauge-invariant. While naively the fact that a = —1 for Maxwell fields and a = for 
gravitons would seem to qualify these theories as exceptional, we expect that as in flat space 
gauge-invariance will prevent the self-energy from acquiring singularities at these values. A 
discussion of such issues will appear soon [72] . It would be very interesting to understand if 
our definition could be further extended to perturbative string theory 16 . 

As a final comment, one would also like to understand any relation of our S-matrix to 
quantities of interest in dS/CFT [73, 74, 75]. Since our S-matrix focuses on asymptotic 

15 This is expected for the Maxwell field but is a matter of significant controversy for the graviton. See 
e.g. [68, 69, 70, 71]. 

16 We thank Eva Silverstein for raising this issue. 
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quantities near future and past infinity, one may expect some relation to exist. However, 
because the Hartle-Hawking state is not itself dual to a particular CFT with fixed sources, 
it is in some sense a derived object in dS/CFT as opposed to a fundamental one [75, 76]. 
Our reliance on the Hartle-Hawking state will thus make any connection somewhat indirect. 
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A KG mode linearization 



This appendix gives some details of the function p ai a 2 L 1 L 2 (l Jj i K) appearing in the lineariza- 
tion formula (2.24). As noted in the main text, the results below may be derived using the 
methods of appendix A of [8], and in particular techniques for manipulating Mellin-Barnes 
representations of hypergeometric functions as described in [77, 41]. 
We introduce a := (D — l)/2 as well as the notation 



ai, a 2 , . . . , ak 
&i,6 2j ...,bj 



r(ai)r(a 2 ) . . .r(a fc ) 



L — cr + z/, L + a + 2a + L — a — z/, —v 

L4 



Kr cos7r(cr + a) C ^ VT T ^ ^-1/2 

H ^ : = 2 L +a7T {T[L-a,L + a + 2a]} 1/2 

and finally the terminating, "regulated" hypergeometric series 

a, 6, -p 

c, d 



3F2 



1 



{o)n{t>)n(-p)r 



71=0 



[1 + n, c + n, d + n] ' 



(A.l) 
(A-2) 
(A-3) 



(A.4) 



with (a) n the Pochhammer symbol. This series is an entire function of each parameter 
a, 6, c, d. 
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The linearization kernel may be then expressed as the double Mellin-Barnes integral 



Pg\g 2 L\L 2 



(», K) = (-lp+ i2 -*)/V +Ll+L W CTlLl A/; 2L2 {r [K - /x, K + 2a + /x]} 



1/2 



/ / S ^<riLi(fl - ^l)^ 2 L 2 (^2 - ^2) 

J V\ J V 2 I 



xr 



X 3^2 



/X - I/i - ZAj, 1 - fJL + Vl + V 2 , 1+g ~ 2 Ll ~ L2 + a + V X + ZAj 
- I* - 1* 

i+A:-L 1 -L 2 +a + yi + y2) l _fE±LrH^ +Vl + v ^ 
1 + i ^~ L 9 1 ~ Z/2 + 2a + // + 1/1 + i/ 2 , 1 + ^" L r Lz - // + fi + ^ 



K-L x -L 2 



; 1 



(A.5) 



The well-known asymptotic behavior of the Gamma function may be used to show that the 
Mellin-Barnes integrals are absolutely convergent. Equation (2.26) follows from the identity 

We are primarily interested locating singular points of p ai a 2 L 1 L 2 (l Jj ^) as a function of 
ji. The singularities of the double MB integral in (A.5) may be determined using standard 
techniques presented, e.g., in Ch. 4 of [78]. Within the integrand there are potentially poles 
due to the Gamma functions r(/i — v\ — v 2 ) an d T(l — /1 + v\ + z/ 2 ). However, the poles in the 
latter Gamma function are canceled by zeros in each term of the hypergeometric series, so 
this Gamma function does not in fact yield singularities. Thus the only source of singularities 
from the double integral is the Gamma function T(/i — v\ — The Vi integrals have poles 
at Vi — di — n and Vi — — [a + 2a) — n, n G No, respectively, so it follows that the double 
integral contributes poles to /1 at poles at precisely the locations (2.25). 

The observant reader will note that the first line of (A.5) contains combinations of gamma 
functions which, when combined with an identical factor in /^(ry), contributes additional 
poles in the integrand at 



ji = K + n, p = —2a — K — n. n G Nq. 



(A.6) 



We did not mention these poles in section 2.2 for two reasons: first, they lie outside the 
strip in the complex \i plane in which we utilize p ai a 2 L 1 L 2 (l Jj ^ K)\ second, these poles do not 
contribute to the behavior of (2.24) at large \r]\ 3> 1, as we now show. Recall the asymptotic 
form (2.14) of f a L(v) which we re- write here as 



U(\v\ > 1) = K 



oh 



L — a, a — a, 2(a + a) 

[l + 0(r]- 2 )], \r]\^>\L-a\. 



(A.7) 
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When \rj\ 3> L\ + L 2 , inserting (A. 7) into (2.24) yields 
fa^J^^dvl > L x +L 2 ) = / 2(fi + a)p (J1(J2LlL2 (fi,K)^ 



K — /x, a — ji ) 2(/i + a) 



r 



+ + 2«, — h « + /i, — 2(/i + «) 

Zi 



-(p+2a) 



[1 + 0(?7" 



(A.8) 



Notice that the factor M^k contains the precisely the combination of gamma functions 
needed to cancel the similar factor in p ai a 2 L 1 L 2 (l Jj i K)- Upon inspection one finds that 
Pa 1 a 2 L 1 L 2 (l Jj i K)M^k contains zeros which cancel the poles due to T(±2(/x + a)). To ob- 
tain the asymptotic form of f ai L 1 (v)fa 2 L 2 (v) from (A.8) we close the integration contour 
with an arc at infinity in the left-half plane for the first term, and with an arc at infinity in 
the right-half plane for the second term, and use Cauchy's integral formula to relate these 
two contour integrals to the sum of residues obtained. The second term in (A.8) yields zero 
because the integrand contains no poles in the right-half plane. The first term in (A.8) 
acquires residues from the poles (2.25), and these poles precisely reproduce the asymptotic 
expansion of f^L^fa^iv)- 



B The explicit form of good operators 

This appendix explores a rather brute force but very explicit method of constructing good 
operators in the context of the model theory (6.1)-(6.41). The basic idea is that, at each 
order, one simply chooses a perturbative field redefinition that removes the unwanted IR- 
divergent contributions but which leaves the behavior at the desired pole unchanged. We 
find such explicit results enlightening, though they become cumbersome at higher orders. 
As a result, we work here only to O(g). The results agree completely with those of section 
7. 

Recall that the difficulty with the naive prescription (3.5) was that the rjo — >> oc limit of 

J S 77=770 

is well-defined only when the correlator satisfies 

(<fo (x) . . .) = (homogeneous) + O^"^ 3 ^" 1 )), rj > 1. (B.2) 

Here (homogeneous) denotes solutions to the homogeneous Klein-Gordon equation (□ — 
M|)/(x) = 0. We refer to the rest of (j>s(x) as the inhomogeneous part. It is straightforward 
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to perturbatively compute the asymptotics of the inhomogeneous part from the Schwinger- 
Dyson equations. For example, at O(g) we have 



(□3-M3 2 ) {<k{x3)<hfa)Mxi)) 



(i) _ 



(0201 (X3) 02 (X 2 ) 01 (Xi ) ) (0) 



= -gW 2 (x 3 ,x 2 )W 1 (x 3 ,x 1 ) 

0(v 3 2+ai ) + 0{rf z 2 - ai - {D - l) ) + 0(7] 3 " <T2+<T1 " (D " 1) ) 

+ 0(t 73 --- ct1 - 2(d - 1) )] [1 + 0( % - 2 )] for \ m \ » 1 (B.3) 

where the final equality summarizes the asymptotic behavior of the Wight man functions. 
Because the action of (D 3 — M|) does not change the asymptotic behavior of power laws of 77, 
the right hand side of (B.3) also describes the inhomogeneous part of (03(x3)02(#2)0i(#i)) - 
This violates the condition (B.2) for —(D — l) < Re (01 + a 2 + 03), and moreover contributes 
unwanted information to (B.l). 

Acting with the R a projector of section 3.3 is thus equivalent (in the sense of section 3.3) 
to constructing some $s(x) whose asymptotics contain the same homogeneous part as (j)$(x) 
but which lacks the divergent inhomogeneous parts. So for —(D — l) < Re (01 + a 2 + 03) 
let P 3 = {(J2 + 0"i, . . . } be the set of weights belonging to terms in (B.3) violating (B.2) and 
consider 



(x)+g 



(B.4) 



where G v {x) denotes an operator with asymptotic behavior 0{rf) + 0(rj ( (73+D x )) at 0(g°). 
We choose the G v {x) to be normalized so that 



M 2 (p) 



Ml 



'201 (a;] 



+ 0( V 



-(r7 3 + D-l)> 



(B.5) 



where (j>{x)\ p denotes the 0(r] p ) part of cj>{x) at O(g ) and M 2 (p) is the squared-mass which 
gives rise to this fall off at 0(g°). The choice of operators G v {x) is not unique, but a natural 
choice is 



7 (x) 



M 2 (p)-Mi 



n 

.qeP 3 ,q^P 



(M 2 (p) - M 2 (q)) 



0201 (x) 



6 2 1 (x) +0{t 1 -^ s+d -V). 
\p 



M 2 (p) - Ml 

With this choice for the & v {x) the new operator $%{x) may also be written 

$3(2) = cf) 3 (x)+gY^ c k n k (f) 2 (f> l (x), 



(B.6) 



(B.7) 



k=0 



with appropriate coefficients c/~. A similar procedure should be applied to $2 0*0 and $i(x). 

It is sufficient to simply change variables from (j>i(x) to <&i(x) in the generating functional 
(the path integral). This is most easily accomplished with <&i(x) defined as in (B.7). The 
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lagrangian for the <&i(x) contains a quadratic part in which each <$>i(x) has a canonically 
normalized kinetic term and a mass term that agrees with that of the original field (j)i{x). 
At order 0{g) the lagrangian contains the cubic interaction c^$3$ 2 $i(^), now with an ad- 
ditional coefficient c, as well as higher-derivative cubic interactions. These higher-derivative 
interactions, taken on their own, produce correlators whose asymptotics have similar fall-off 
properties to the simple cubic interaction. But these higher-derivative interactions conspire 
to cancel the unwanted inhomogeneous terms in (B.3) without altering the homogeneous 
part. One could say that the new higher-derivative interactions are "IR counterterms," 
though we stress that the theory has not been altered nor does the theory need such terms 
in order for correlation functions to be well-defined. 

We can now compute 0(g) amplitudes. To keep the equations tractable let us consider 
the case where there is only one problematic term at 0(7y cri+cr2+(73 ). In this case the <&i(x) 
defined in (B.6) are 



$ 3 Or) 



(x) + 



9 



M 2 (a 2 + (r 1 )-Ml 



020lO) 



= 03 0) + 



"0201 (x) 



and likewise for $2 0*0, $ 



Ax 



(cr 3 - Gi - cr 2 )(cri + a 2 + cr 3 + D - 1) 
The 0{g) 3-pt. function of the corrected operators is 



(B.8) 



{^{X^ 2 (x 2 )^ 1 (x l )f ) = (0 3 (X 3 )02(X 2 )01(X1)) (1) 

9 



M\a 2 + <n) 


-Ml 


+ 9 






-Ml 


+ 9 






-Ml 



(0201 (^3)02(^2)01 (Xi)) (0) 
(03 (^3)0301 (X 2 )<t>l (^l)) (0) 
(03(£3)02(£2)0302(£l)) (O) - 



(B.9) 



The 1 — > 2 transition amplitude is then 

/ (n 3 n 2 |ni)f ) = ig ul^u^y) 
ig I fiE !/ 1 



U 2 U\ 



-ig \ <fi? 



2 (M 2 (a 2 + <7i) - Ml) 
u* 3 u* 2 (x) 



^ u u* 3 (x) + (2 o 3) 



V=+Vo 



Ui 



M 2 (a 3 + a 2 )-Ml 



(B.IO) 



-770 



Here we have regulated right-hand-side by imposing \rj\ < rjo since each term is individually 
divergent as 1 770 1 00. We have also explicitly symmetrized the operators of the final state. 

The first term is just the regulated version igIi(r]o) of the integral considered in §6.1. 
Recall that our computation of h(r]o) for heavy fields used the linearization formula (2.22) 
on UsU2(x) then deformed the ji integration contour to Re fi < — Re(ai + D — 1). In the 
process we acquired residues due to simple poles at \i = G\ and (ai + D — 1). A similar 
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procedure can be used in the present case but now, since — Re(<7i + D — 1) < Re(<72 + 03), 
deforming the contour to the left picks up an additional residue from the simple pole in the 
linearization kernel at fi = a 2 + 03. This new pole now provides the leading behavior of 
u 3 u 2 (x). Thus 



7i(t7o) = (RHS of eq. 6.19) 
1 



M 2 (02 + 0-3)- Ml 
+ ..., 



dE"(y) 



* * 

u 3 u 2 



ri=Vo 



u 3 u 2 



V^iii(y) 



V=-Vo 

(B.ll) 



where the ellipses denote the remaining integral over some contour with Re \i < — Re(o"i + 
D — 1). Contributions from this remaining contour integral vanish in the limit 770 —> 00. 

We now examine the terms on the second and third lines of (B.10). Note that only the 
leading term of each KG mode Ui{x) = 0(r] ai ) will contribute to these surface integrals in 
the limit rj — >> 00. Noting that 



1 



(cr 3 - (J 2 - (Ti) 



Ml 



(2 <-> 3) 



2 LM 2 (<72 + ^l) 

it follows that the contribution of the integrals at 77 



(03 + 02 ~ 0l) 

M 2 (0 3 + 2 ) -M 2 ' 
+H may also be written as 



(B-12) 



ig / dE 1 ' 



U 2 U\ 



2(M 2 (0 2 + 0i)-M 2 ) 



V,^(x) + (2 o 3) 



'?=+% 



M 2 (0 3 + 2 ) - Ml 



v u\u* 2 {y) 



+ ... 



V=+Vo 



(B.13) 



where again ellipses denotes terms that vanish as rjo oc. It follows that for t]q —> oc the 
last two lines of (B.10) precisely cancel the new residue term in h(r]o). Thus the amplitude 
continues to be given by ig times (6.19). 

It will not surprise the reader that the same result turns out to be valid for all > 0. A 
detailed check shows that similar cancellations occur in the presence of further problematic 
terms. For every divergent term 0(r] (Jl+(72+(73 ), 0(r] (7l+(72+(73 - 2 ), O^-^i+^+^-p-i)^ etc ^ 
one acquires a set of surface integrals arising from the residue of a pole in the derivation of 
h(vo)- These surface integrals are exactly canceled by the contribution of the appropriate 
G v {x) operators. It is also straightforward (if tedious) to verify similar cancellations of the 
divergent terms in the remaining 0(g) amplitudes. In agreement with section 7, we see that 
the results are again given by the corresponding final formulae in section 6.3. 
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